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FOREWORD 


The  efficient  application  of  turbojet  engines  to  the  propulsion  of  aircraft 
at  high  supersonic  Mach  numbers  requires  extensive  knowledge  of  the  characteristics 
of  inlets  and  induction  systems.  In  this  connection,  the  Research  Deportment  of 
United  Aircraft  Corporation  has  been  investigating  various  turbojet  inlet  design 
parameters  under  Contract  NOa(s)  55-133-c  with  the  Power  Plant  Division  of  the 
Bureau  of  Aeronautics,  Department  of  the  Navy,  to  provide  basic  Information  for  the 
design  of  high  Mach  number  inlets.  The  results  of  the  investigations  performed 
under  Lots  I  and  II  of  this  contract  are  summarized  in  Refs.  1  through  6.  Additional 
investigations  performed  under  Lot  III  of  the  same  contract  pertain  to  the  following 
subjects: 


1.  Performance  of  a  Mach  1+  .0  variable -geometry  axi  symmetric  inlet  having 
external -plus -internal  compression. 

2.  Efficiency  of  flush  oblique  nozzles  exhausting  into  supersonic  streams 
having  Mach  numbers  up  to  4.0. 

3.  Turbulent  boundary  layer  characteristics  in  supersonic  streams  having 
adverse  pressure  gradients. 


This  report  summarizes  the  results  of  the  studies  of  the  effect  of  adverse 
pressure  gradients  on  boundary  layer  characteristics.  The  investigations  which 
pertain  to  the  remaining  subjects  are  summarized  in  Refs.  7  and  8. 
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Turbulent  Boundary  Layer  Characteristics  in  Supersonic  Streams 
Having  Adverse  Pressure  Gradients 

SUMMARY 


Investigations  were  conducted  to  determine  the  thickness  and  profile  shape 
characteristics  of  turbulent  boundary  layers  in  supersonic  streams  having  adverse 
pressure  gradients  similar  to  those  which  are  encountered  in  supersonic  inlets.  The  ' 
program  included  tests  conducted  at  Mach  numbers  from  2.0  to  3-5  using  two-dimensional 
and  axisymmetric  curved- surface  models  to  produce  the  adverse  pressure  gradients. 

The  magnitude  of  the  gradients  relative  to  the  boundary  layer  thickness  at  the  be¬ 
ginning  of  the  gradient  was  varied  by  employing  models  having  different  radii  of 
curvature  and  by  changing  the  boundary  layer  thickness  at  the  beginning  of  the 
gradient.  The  over-all  pressure  rise  in  most  cases  was  greater  than  the  value  which 
would  cause  a  turbulent  boundary  lay~r  to  separate  if  the  pressure  rise  were  created 
by  an  oblique  shock  wave.  Some  tests  were  conducted  to  determine  the  effect  of 
boundary  layer  removal  from  the  curved  surfaces  on  boundary  layer  growth  and  on  the 
pressure  gradients  which  could  be  realized  without  separation.  An  analytical  in¬ 
vestigation  was  also  conducted  so  that  the  results  of  the  experimental  Investigation 
could  be  applied  to  the  prediction  of  cases  outside  the  range  of  the  experiments. 


1  the 


It  is  shown  that  the  boundary  layer  momentum  thickness  can  be  predicted  from 
von  Kamian  boundary  layer  momentum  equation  kbrrb-fthat  the  measured  values  of  boundary 
layer  profile  shape  arc  ij^poor  agreement  wlth^tttJb  values  computed  from  procedures 
derived  by  extendinft^rTWent ional—il^ l  -predicting  profile  shape  in  subsonic 

flow.  It  was  found'OfaaO  boundary  layers  in  adverse  pressure  gradients  in  supersonic 
streams  are  much  less  susceptible  to  separation  than  is  Indicated  by  extension  of 
conventional  procedures  developed  for  subsonic  flow;  therefore,  supersonic  inlets  can 
be  made  much  shorter  and  lighter  than  inlets  which  are  designed  on  the  basis  of  these 
conventional  procedures.  A  new  procedure  for  calculating  boundary  layer  profile 
shapes,  developed  in  this  report,  is  shown  to  provide  a  good  correlation  between  ex¬ 
perimental  and  calculated  values  of  boundary  layer  profile  shapes  in  adverse  pressure 
gradients  created  by  curved  surfaces.  This  new  procedure  is  based  on  the  experimental 
observation  that  the  station  at  which  high-energy  free- stream  flow  actually  mixes  into 
a  turbulent  boundary  layer  in  an  adverse  pressure  gradient  is  well  downstream  of  the 
station  at  which  flow  would  have  to  mix  in  order  to  maintain  a  flat-plate  profile. 


This  project  was  undertaken  for  the  Department  of  the  Navy,  Bureau  of  Aeronautics 
under  Lot  III  of  Amendment  No.  7  to  Contract  N0a(s)  55-133-c. 
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CONCLUSIONS 


1.  Boundary  layers  in  adverse  pressure  gradients  in  supersonic  streams  are 
much  less  susceptible  to  separation  than  is  indicated  by  extension  of  previously 
developed  procedures  for  subsonic  flowj  therefore,  supersonic  inlets  can  bejnade 
much  shorter  and  lighter  than  inlets  which  are  designed  on  the  basis  of  extrapolated 
subsonic  theory. 

2.  The  shape  of  a  turbulent  boundary  layer  in  adverse  pressure  gradients 
generated  by  curved  surfaces  in  supersonic  flow  can  be  predicted  with  much  greater 
accuracy  by  the  lag-length  procedure  developed  in  this  report  than  by  existing 
procedures  based  on  extrapolation  of  data  from  subsonic  flow.  The  lag-length  pro¬ 
cedure  makes  use  of  experimental  data  which  show  that  high-energy  free- stream  flow 
begins  to  mix  with  the  low-energy  boundary  layer  flow  at  a  station  located  at  a 
distance  (called  the  lag  length)  from  0  to  60  boundary  layer  momentum  thicknesses 
downstream  of  the  station  where  flow  would  have  to  mix  to  maintain  a  flat-plate 
profile . 


3.  Dimensionless  lag  lengths  obtained  by  dividing  the  measured  lag  lengths  by 
the  boundary  layer  momentum  thickness  are  a  function  of  Reynolds  number  but  are  in¬ 
dependent  of  Mach  number  and  pressure  gradient  for  a  wide  range  of  Mach  numbers, 
pressure  gradients,  and  Reynolds  numbers  for  curved  surfaces  in  supersonic  flow. 

4.  The  momentum  thickness  of  a  turbulent  boundary  layer  in  adverse  pressure 
gradients  generated  by  curved  surfaces  in  supersonic  flow  can  be  calculated  from  the 
von  Karman  boundary  layer  momentum  equation. 

5-  The  use  of  boundary  layer  bleed  results  in  a  large  increase  in  the  pressure 
gradient  which  can  be  tolerated  by  a  turbulent  boundary  layer  in  supersonic  flow 
without  separating. 

6.  The  decrease  in  the  momentum  thickness  of  a  turbulent  boundary  layer  due  to 
bleed  is  approximately  half  that  calculated  with  the  assumption  that  the  bleed  flow 
has  zero  axial  velocity. 
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INTRODUCTION 


The  characteristics  of  boundary  layers  in  adverse  pressure  gradients  at  super¬ 
sonic  speeds  have  a  particularly  large  influence  on  the  performance  of  supersonic  in¬ 
lets.  The  thickness  of  the  boundary  layer  approaching  the  normal  shock  at  the  throat 
has  a  direct  influence  on  the  pressure  recovery  through  the  normal  shock  (see  Ref.  l). 
In  addition,  boundary  layer  separation  can  prevent  attainment  of  a  low  Mach  number 
upstream  of  the  normal  shock  and,  hence,  prevent  attainment  of  high  inlet  pressure 
recoveries.  If  the  inlet  is  made  longer  than  the  minimum  length  necessary  to  prevent 
separation,  unnecessary  structural  weight  will  result. 

The  object  of  the  research  described  herein  was  to  establish  criteria  for  pre¬ 
dicting  the  characteristics  of  boundary  layers  in  supersonic  streams  having  adverse 
pressure  gradients  to  provide  a  more  rational  basis  for  the  design  of  supersonic 
inlets. 


RESUME  OF  PREVIOUS  WORK 


Most  of  the  data  available  in  the  literature  (see,  for  example.  Refs-  9  through 
19)  were  obtained  from  tests  in  which  the  pressure  rise  was  created  by  oblique  or 
normal  shock  waves  (severe  adverse  pressure  gradients)  and  in  which  the  mixing  of 
flow  from  the  free  stream  into  the  boundary  layer  had  little  effect  on  the  shape  of 
the  boundary  layer  profile  within  the  severe  adverse  pressure  gradient.  The  portion 
of  the  data  from  these  reports  dealing  with  the  maximum  pressure  rise  that  can  be 
obtained  without  separation  has  been  explained  in  Refs.  20  through  24  on  the  basis 
that  there  is  little  or  no  mixing  of  free-stream  flow  with  the  flow  in  the  boundary 
layer. 

The  only  data  which  Is  known  to  be  available  on  boundary  layer  characteristics 
in  supersonic  streams  for  adverse  pressure  gradients  more  moderate  than  those  created 
by  shock  waves  is  contained  in  Ref.  25.  This  NASA  memorandum  provides  separation 
criteria  only  and  does  not  include  information  on  the  actual  shape  of  the  boundary 
layer  within  the  adverse  pressure  gradients.  Attempts  to  explain  the  pressure-rise 
data  presented  therein  on  the  basis  of  the  theories  of  Refs.  20  through  24  were  not 
successful  because  the  theories  of  these  five  reports  neglect  the  effect  of  mixing 
the  free-stream  flow  with  the  flow  within  the  boundary  layer.  The  only  known  theory 
for  the  calculation  of  boundary  layer  shape  in  a  moderate  adverse  pressure  gradient 
at  supersonic  speeds,  given  in  Ref.  26,  is  an  extension  of  the  theory  developed  in 
Ref.  27  for  subsonic  flow.  The  latter  reference  presents  one  of  many  similar  theories 
(see  Refs.  28,  29,  and  30)  which  have  been  developed  for  the  prediction  of  boundary 
layer  shape  for  subsonic  flow.  It  is  shown  in  the  present  report  that  the  terms  In 
the  equations  of  Ref.  26  which  account  for  the  mixing  of  free-stream  flow  with  flow 
in  the  boundary  layer  severely  underestimate  the  improvement  in  boundary  layer  shape 
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due  to  mixing  and,  hence,  severely  underestimate  the  pressure  rise  which  can  be  ob¬ 
tained  without  separation  in  moderate  adverse  pressure  gradients. 

The  effect  of  pressure  gradient  on  boundary  layer  thickness  must  be  evaluated 
in  supersonic  inlets  because  the  total  pressure  recovery  through  the  normal  shock  at 
the  throat  of  the  inlet  is  a  function  of  the  thickness  of  the  boundary  layer  approach¬ 
ing  the  throat  (see  Ref.  l).  Although  some  data  are  available  on  boundary  layer 
growth  in  severe  adverse  pressure  gradients  (see,  for  example,  Refs.  9,  l4,  and  19), 
no  data  are  known  to  be  available  on  boundary  layer  growth  in  moderate  adverse  pressure 
gradients.  In  calculating  the  thickness  of  a  boundary  layer  in  an  adverse  pressure 
gradient,  it  is  usually  convenient  to  employ  the  boundary  layer  momentum  thickness. 

Most  procedures  for  calculating  boundary  layer  momentum  thickness  in  adverse  pressure 
gradients  are  based  on  the  use  of  the  von  K^rm4n  boundary  layer  momentum  equation. 

Two  adaptations  of  the  von  K^trmAn  boundary  layer  momentum  equation  for  calculating 
momentum  thickness  growth  in  an  adverse  pressure  gradient  at  supersonic  speeds  are 
given  in  Refs .  26  and  31 • 


EXPERIMENTAL  INVESTIGATION 


Description  of  Test  Equipment 


Tunnel 


Tests  to  determine  boundary  layer  characteristics  in  adverse  pressure  gradients 
were  conducted  in  the  UAC  Research  Department  4.5  x  4. 5- in.  continuous-flow  super¬ 
sonic  wind  tunnel.  Dry  air  was  supplied  to  the  tunnel  at  a  total  pressure  of  approxi¬ 
mately  one  atmosphere  and  a  total  temperature  of  approximately  150  F.  The  Reynolds 
number  per  in.  in  the  tunnel  test  section  varied  from  0.27  x  10^  at  a  Mach  number  of 
2.0  to  0.12  x  10^  at  a  Mach  number  of  3*5*  Two  test-section  configurations  were 
employed,  one  for  tests  of  two-dimensional  models  and  one  for  tests  of  axisymmetric 
models  (a  description  of  the  models  is  given  in  the  following  subsection). 

A  sketch  and  a  photograph  of  t1  e  test-section  configuration  used  in  tests  of  the 
two-dimensional  curved-surface  model  are  presented  in  Figs.  1  and  2,  respectively. 

A  flat  surface  on  the  centerline  of  tie  tunnel  (centerline  block)  produced  the 
boundary  layer  approaching  the  models.  The  use  of  this  centerline  block  reduced  the 
height  of  the  test  section  by  one  half  to  the  dimensions  given  in  Fig.  1.  The  use 
of  surface  roughness  near  the  tunnel  throat,  and  the  injection  of  air  into  the 
boundary  layer  through  the  injection  holes  (see  Fig.  l)  were  employed  to  trip  the 
boundary  layer  to  insure  turbulent  flow.  The  flow- injection  system  was  also 
utilized  as  a  method  of  increasing  the  boundary  layer  thickness  at  the  upstream  end 
of  the  curved  surface.  The  flow  rate  into  the  plenum  and,  hence,  the  thickness  of 
the  boundary  layer  approaching  the  two-dimensional  models  was  regulated  by  a  valve 
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vented  to  atmosphere. 

Provision  was  made  for  boundary  layer  bleed  from  the  surface  of  selected  two- 
dimensional  curved- surface  models  into  the  bleed  plenum  shown  in  Fig.  1.  The  flow 
passing  through  the  bleed  plenum  was  ducted  through  a  flow-measurement  bellmouth  and 
a  control  valve  to  a  region  of  low  static  pressure  in  the  tunnel  downstream  of  the 
models . 

The  test-section  conf iguration  employed  in  tests  of  the  axisymmetric  curved- 
surface  models  was  obtained  by  replacing  the  centerline  block  with  a  nozzle  block 
which  matched  the  top  nozzle  block  shown  in  Figs.  1  and  2.  The  axisymmetric  models 
were  mounted  in  an  adjustable  support  which  was  attached  to  the  tunnel  side  wall. 

The  probe  employed  to  measure  the  boundary  layers  (see  Figs.  1  and  2)  was  driven 
by  a  mechanism  which  allowed  the  probe  to  traverse  at  any  angle  in  the  vertical  plane 
so  that  the  traverse  could  be  made  normal  to  the  model  surface.  The  position  of  the 
probe  and  the  magnitude  of  the  probe  pitot  pressure  were  measured  by  electrical 
transducers  and  the  output  from  these  transducers  was  recorded  on  IBM  punch  cards 
at  selected  values  of  probe  position. 

Models 


Three  types  of  two-dimensional  curved- surface  models  which  did  not  employ 
boundary  layer  bleed  were  utilized:  one  type  was  a  reference  flat  plate,  the  second 
type  had  a  constant  radius  of  curvature  over  the  entire  curved  surface  (see  Figs.  3 
and  4),  and  the  third  type  had  two  different  radii  of  curvature  (see  Figs.  5  and  6). 
The  two-radii  curved- surface  models  were  formed  by  the  addition  of  extensions  to  the 
single-radius  curved- surface  model  having  a  radius  of  curvature  of  6  in.  and  total 
turning  angle  of  20  deg.  Extensions  were  available  which  had  radii  of  1.5  and  3  in. 
and  which  increased  the  total  turning'  angle  to  30  or  35  deg. 

Provisions  were  made  for  boundary  layer  bleed  through  perforations  on  selected 
single-radius  two-dimensional  curved-surface  models.  These  models  are  shown  in 
Figs.  3  and  4.  The  perforations,  which  were  l/l6  in.  in  diameter  and  located  on 
l/4-in.  centers  across  the  model,  were  distributed  over  the  entire  curved  surface  of 
the  model.  The  porosity  of  the  model  could  be  varied  by  filling  the  perforations 
with  plastic. 

The  sides  of  all  models  having  total  turning  angles  greater  than  20  deg  were 
cut  back  from  the  walls  of  the  tunnel  in  order  to  minimize  the  effect  of  tunnel  side- 
wall  boundary  layer  on  the  flow  over  the  models.  The  models  with  modifications  to 
provide  this  change  in  geometry  are  shown  in  Figs.  3,  4,  5,  and  6.  A  bypass  angle, X, 
of  5  deg  was  employed  for  all  models  having  total  turning  angles  greater  than  20  deg. 
However,  in  one  test  with  a  model  having  a  total  turning  angle  of  20  deg,  a  check 
on  model  performance  was  made  with  bypass  angles  of  0  and  20  deg. 
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A  sketch  and  a  photograph  of  the  axisymraetric  models  are  given  in  Figs.  7  and  8. 
Single-radius  models  having  radii  of  curvature  of  1-5  and  3  in.  and  total  turning 
angles  of  20  and  30  deg  were  employed.  The  grooves  and  grit  shown  in  Figs.  7  and  8 
were  employed  to  promote  boundary  layer  transition  and  control  the  boundary  layer 
thickness  at  the  beginning  of  the  curved  surface  for  most  of  the  tests.  For  one 
test,  however,  a  perforated  section  was  employed  upstream  of  the  curved  surface  (see 
photograph.  Fig.  8)  to  increase  the  boundary  layer  thickness  by  bleeding  air  from 
inside  the  model  into  the  boundary  layer. 

The  model  configurations  and  test  conditions  for  which  data  are  reported  are 
listed  in  Table  I. 


Data  Reduction  Procedure 

Conventional  data  reduction  procedures  employed  in  the  calculation  of  boundary 
layer  parameters  on  flat  plates  cannot  be  used  for  calculating  parameters  on  curved 
surfaces  because  large  static  pressure  gradients  normal  to  the  curved  surfaces  are 
required  to  provide  the  centripetal  force  necessary  to  turn  the  flow.  Examples  of 
the  measured  and  calculated  pressure  distributions  in  a  boundary  layer  having  a 
large  static  pressure  gradient  normal  to  the  wall  are  given  in  Fig.  9-  Although  the 
wall  static  pressure  was  measured,  the  static  pressure  distribution  through  the 
boundary  layer  could  not  be  determined  accurately  by  use  of  a  static  pressure  probe 
because  of  local  flow  angularity.  Therefore,  it  was  necessary  to  calculate  the  local 
static  pressure  distribution  using  measured  values  of  pitot  pressure,  wall  static 
pressure,  and  free-stream  total  pressure.  The  first  approximation  to  the  variation 
of  local  static  pressure  (Pgf)  with  distance  from  the  wall  was  obtained  from  the 
measured  local  pitot  pressures  and  a  total  pressure  equal  to  free-stream  total 
pressure.  An  example  of  this  variation  of  Pg^,  (divided  by  Pto)  is  shown  in  Fig.  9* 
The  values  shown  are  equal  to  the  true  local  static  pressure  outside  the  boundary 
layer  where  the  local  total  pressure  is  equal  to  free-stream  total  pressure,  but  are 
incorrect  within  the  boundary  layer  where  the  local  total  pressures  are  lower  than 
free-stream  total  pressure.  In  the  next  step  of  the  calculation  procedure,  it  was 
assumed  that  the  values  of  Fg^.  were  correct  at  a  value  of  y  equal  to  yMATCH  where 
yMATCH  was  determined  from  an  arbitrary  constant,  K,  times  the  value  of  y  at  the 
maximum  value  of  pitot  pressure.  The  variation  of  static  pressure  from  y  =  0  to 
y  =  yMATCH  was  obtained  by  fitting  a  quartic  curve  between  the  wall  static  pressure 
and  the  value  of  Ps^.  at  y  =  yMATCH-  Since  a  quartic  curve  is  determined  by  five 
conditions,  three  additional  conditions  over  and  above  the  static  pressure  at  the 
wall  and  the  value  of  Pg^,  at  y  =  yMATCH  were  needed  to  define  it.  These  three  con¬ 
ditions  were  determined  by  the  requirement  that  the  static  pressure  gradient  at  the 
wall  be  zero  (determined  on  the  basis  of  a  velocity  of  zero  at  the  wall)  and  that 
the  first  and  second  derivatives  of  the  static  pressure  at  y  =  yMATCH  equal  to 
the  values  determined  from  the  Pg^.  curve. 
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The  local  total  pressure  in  the  boundary  layer  can  be  calculated  from  the  local 
static  pressure  and  the  measured  pitot  pressure.  An  example  of  the  variation  of  the 
local  total  pressure  in  the  boundary  layer  with  distance  normal  to  the  surface  is 
given  in  the  upper  part  of  Fig.  9.  The  boundary  layer  thickness,  8 ,  is  defined  as 
the  value  of  y  at  which  the  local  total  pressure  is  equal  to  95$  of  the  tunnel  total 
pressure.  For  flow  with  no  static  pressure  gradient  normal  to  the  wall  and  a  Mach 
number  outside  the  boundary  layer  of  2.5,  this  definition  of  the  boundary  layer 
thickness  resulted  in  a  velocity  at  the  outside  edge  of  the  boundary  layer  which 
was  99. of  the  velocity  outside  the  boundary  layer. 

It  is  conventional  in  the  calculation  of  boundary  layer  parameters  on  a  flat 
plate,  where  the  conditions  outside  the  boundary  layer  are  uniform,  to  base  the  dis¬ 
placement  thickness  and  momentum  thickness  of  the  boundary  layer  on  the  conditions 
outside  the  boundary  layer.  However,  in  the  tests  conducted  in  this  program,  the 
variation  of  conditions  outside  the  boundary  layer  with  distance  from  the  wall  was 
so  large  that  this  method  no  longer  provided  a  convenient  set  of  reference  conditions. 
Therefore,  it  was  decided  to  define  boundary  layer  displacement  thickness  and  boundary 
layer  momentum  thickness  on  the  basis  of  conditions  calculated  from  the  wall  static 
pressure.  The  areas  represented  by  the  integrals  Involved  in  the  definition  of  the 
boundary  layer  quantities  are  shown  in  Fig.  10. 

A  value  of  match  ratio,  K  =  1.3,  wa3  chosen  for  reduction  of  the  data.  Calcu¬ 
lations  showed  that  an  increase  of  K  from  1.2  to  1-5  (a  change  in  yfoATCH  of  25$) 
resulted  in  an  increase  in  boundary  layer  thickness,  8 ,  of  approximately  5$  and  a 
decrease  in  boundary  layer  displacement  thickness,  8*,  and  momentum  thickness,  9 , 
of  only  3$.  Therefore,  the  error  introduced  by  the  arbitrary  choice  of  K  was  not 
significant  in  the  evaluation  of  8,  8*,  or  9 . 


Reference  Flat  Plate  Studies 

The  uniformity  of  the  flow  in  the  test  region  with  the  reference  flat  plate  in¬ 
stalled  in  the  tunnel  is  illustrated  in  Fig.  11.  Although  the  Mach  number  indicated 
by  the  wall  static  pressure  at  the  tangent  points  for  the  different  models  (Station  0) 
varied  from  model  to  model,  the  tunnel  Mach  number  referred  to  in  presenting  data  is 
the  nominal  test-section  Mach  number  (2.0,  2.5,  3.0,  or  3*5).  Typical  values  of 
boundary  layer  characteristics  along  the  centerline  of  the  reference  flat  plate  are 
given  in  Fig.  12. 

The  effect  of  Mach  number  on  the  boundary  layer  shape  parameter,  H,  for  a  1/7- 
power  velocity  profile  is  given  in  Fig.  13.  Because  of  the  large  variation  of  H 
with  Mach  number  for  l/7-power  profiles,  the  absolute  numerical  value  of  H  in  com¬ 
pressible  flow  loses  the  significance  that  it  has  in  incompressible  flow,  where  H 
for  turbulent  flow  on  a  flat  plate  is  near  1.3  and  where  H  at  separation  is  usually 
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between  1.8  and  2.6.  In  order  to  employ  a  boundary  layer  shape  parameter  whose 
absolute  numerical  values  retain  physical  significance  in  compressible  flow,  a 
reduced  boundary  layer  shape  parameter,  G,  was  adopted,  where 


(1) 


Hjj  =  7  in  the  preceding  equation  was  obtained  from  the  tables  in  Ref.  31  (see  Fig.  13). 
Values  of  this  reduced  boundary  layer  shape  parameter  calculated  from  the  data  given 
in  Fig.  12  are  presented  in  Fig.  l4. 

As  noted  in  a  preceding  section,  the  boundary  layer  thickness  approaching  the 
models  was  controlled  by  the  pressure  in  the  injection  plenum  tank  shown  in  Fig.  1. 

The  effects  of  injection  plenum  pressure  on  the  boundary  layer  thicknesses, 8 ,  8*, 
and  6 ,  and  on  the  boundary  layer  shape  parameters,  H  and  G,  are  shown  in  Figs.  15  and 
16 .  There  was  a  sudden  change  in  boundary  layer  characteristics  at  an  injection 
plenum  pressure  equal  to  approximately  5$  of  the  tunnel  total  pressure  because  an 
injection  plenum  pressure  equal  to  or  greater  than  this  value  was  required  to  fix 
boundary  layer  transition  at  a  point  near  the  injection  station.  The  effect  of  air 
injection  on  the  boundary  layer  velocity  profiles  for  the  reference  flat  plate  is 
shown  in  Fig.  17-  Although  air  injection  resulted  in  an  increase  in  boundary  layer 
thickness  by  a  factor  greater  than  two  (see  Fig.  15),  it  resulted  in  very  little 
change  in  the  shape  of  the  boundary  layer  velocity  profiles  in  Fig.  17 . 

The  effect  of  Reynolds  number  on  the  reduced  boundary  layer  shape  parameter,  G, 
for  the  reference  flat  plate  is  given  in  Fig.  18.  The  line  faired  through  these 
data  points  (noted  as  Gp)  is  used  as  a  reference  value  of  reduced  boundary  la;  er 
shape  parameter  for  some  of  the  analyses  which  are  discussed  in  the  following 
sections.  The  qualitative  variation  of  G  with  Reynolds  number  is  expected  on  the 
basis  of  the  variation  of  the  exponential  profile  parameter,  N,  with  Reynolds  number 
as  given  in  Ref.  32.  However,  the  absolute  values  of  G  shown  for  any  given  value 
of  R g  in  Fig.  l8  are  approximately  0.03  greater  than  the  values  of  G  determined  from 
the  tables  in  Ref.  31  using  the  values  of  N  given  in  Ref.  32.  The  reason  for  this 
discrepancy  is  that  the  values  of  N  given  in  Ref.  32  were  determined  from  faired 
lines  which  fit  the  velocities  in  the  outer  portion  of  the  boundary  layer  but  did 
not  fit  the  measured  velocities  near  the  wall. 


Boundary  Layer  Characteristics  Without  Bleed 

Thirty-eight  runs  were  made  with  models  having  no  boundary  layer  bleed,  a  nan 
being  defined  as  the  set  of  measurements  necessary  to  determine  the  variation  of 
boundary  layer  characteristics  with  distance  along  a  curved- surface  model  at  one 
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Mach  number  and  one  value  of  boundary  layer  thickness  approaching  the  model.  Of 
these  38  runs,  21  were  obtained  with  single- radius  two-dimensional  models,  11  with 
two-radii  two-dimensional  models,  and  6  with  axisymmetric  models.  The  values  of 
initial  boundary  layer  momentum  thickness  approaching  the  models  are  given  in 
Table  I.  None  of  these  38  runs  contain  data  for  which  separation  was  indicated  at 
any  point  along  the  curved  surface.  Sample  data  obtained  from  the  testB  are  dis¬ 
cussed  in  the  following  paragraphs . 

Typical  schlieren  photographs  obtained  during  tests  of  single-radius  two- 
dimensional  curved-surface  models  are  presented  in  Fig.  19-  The  curved  shock  waves 
in  the  lower  portion  of  these  schlieren  photographs  were  formed  by  coalescence  of 
the  Mach  waves  originating  from  the  two-dimensional  curved  surfaces.  By  following 
streamlines,  it  was  determined  that  the  total  pressure  loss  due  to  this  oblique  shock 
wave  near  the  outside  edge  of  the  boundary  layer  at  the  downstream  end  of  the  model 
was  less  than  1$  of  free-stream  total  pressure  for  all  models  tested. 

All  models  employed  in  the  tests  discussed  in  this  report  were  designed  so  that 
a  Mach  wave  emanating  from  the  sidewall  of  the  tunnel  at  Station  0  and  progressing 
downstream  at  an  angle  determined  from  the  Mach  number  calculated  from  wall  static 
pressure  would  not  intersect  the  centerline  of  the  model  until  a  point  downstream  of 
the  curved  surface.  This  limitation  on  the  maximum  model  length  did  not,  however, 
eliminate  the  possibility  of  changes  in  the  flow  on  the  centerline  of  the  model  due 
to  propagation  of  end  effects  through  the  boundary  layer.  To  investigate  this 
possibility,  a  series  of  tests  was  conducted  during  the  initial  portion  of  the  ex¬ 
perimental  program  to  determine  the  influence  of  the  boundary  layer  on  the  side  walls 
of  the  tunnel  on  the  flow  over  the  two-dimensional  curved- surface  models  The  effect 
of  wall  boundary  layer  was  altered  by  changing  the  bypass  angle  of  the  model  (see 
Figs.  3  through  6).  During  preliminary  tests  the  models  haying 'total  turning  angles 
of  20  deg  functioned  satisfactorily  with  a  bypass  angle  of  zerd.  However,  for  a 
bypass  angle  of  zero  and  total  turning  angles  greater  than  20  deg,  the  flow  separated 
from  the  models  because  of  the  effect  of  the  adverse  pressure  gradient  on  the  thick 
wind-tunnel  sidewall  boundary  layer.  Although  this  separation  was  eliminated  by  the 
use  of  a  bypass  angle  of  20  deg,  the  static  pressures  near  the  downstream  end  of 
the  models  with  bypass  angles  of  20  deg  were  considerably  less  than  those  predicted 
by  theory,  and  the  static  pressures  on  the  ends  of  the  model  were  considerably  less 
than  on  the  model  centerline.  These  results  indicated  a  severe  sideways  flow  of  air 
within  the  boundary  layer  for  a  bypass  angle  of  20  deg.  However,  it  was  possible 
with  a  bypass  angle  of  5  deg  to  test  models  with  turning  angles  greater  than  20  deg 
without  separation,  with  static  pressure  distributions  approximately  equal  to  those 
predicted  by  theory,  and  with  relatively  flat  distributions  of  static  pressure 
across  the  models.  Therefore,  a  bypass  angle  of  5  deg  was  employed  for  all  subse¬ 
quent  tests  of  models  having  total  turning  angles  greater  than  20  deg. 

Typical  data  showing  the  variation  of  wall  static  pressure  with  distance  across 
a  two-dimensional  curved-surface  model  for  three  different  values  of  initicl 
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boundary  layer  thickness  are  shown  on  the  upper  part  of  Fig.  20.  During  tests  of  one 
model  which  employed  a  5-deg  bypass  angle,  measurements  of  the  boundary  layer 
characteristics  at  several  stations  across  the  model  near  the  downstream  end  of  the 
curved  surface  indicated  that  there  was  no  significant  change  in  boundary  layer  shape 
or  thickness  with  distance  across  the  model. 

Examples  of  the  variation  of  wall  static  pressure  with  distance  from  the  curved- 
surface  tangent  point  of  a  two-dimensional  curved- surface  model  are  given  in  the 
lower  part  of  Fig.  20.  The  small  differences  between  the  measured  and  theoretical 
(based  on  method  of  characteristics)  static  pressure  distributions  shown  in  the  lower 
plot  of  Fig.  20  were  caused  by  the  increased  boundary  layer  growth  as  compared  to  the 
growth  on  the  reference  flat  plate  and  by  gradients  in  the  test  section  (see  Fig.  11 )  ■ 
The  variation  of  Mach  numbers  determined  from  the  wall  static  pressures  given  in 
Fig.  20  with  distance  from  the  curved- surface  tangent  point  is  shown  in  Fig.  21. 

Also  shown  in  Fig.  21  are  Mach  numbers  near  the  outside  edge  of  the  boundary  layer 
determined  from  the  boundary  layer  reduction  procedure  outlined  in  Figs.  9  and  10. 

The  difference  between  the  Mach  number  indicated  by  the  wall  static  pressure  and  the 
Mach  number  at  the  outside  edge  of  the  boundary  layer  is  a  measure  of  the  difference 
in  static  pressure  required  to  provide  the  centripetal  force  necessary  to  turn  the 
flow. 


Typical  boundary  layer  velocity  profiles  for  different  stations  along  a  model 
are  given  in  Fig.  22.  The  velocities  given  in  the  figure  are  referenced  to  the 
nominal  tunnel  velocity  because  of  the  variation  of  velocity  with  distance  from  the 
wall  outside  the  boundary  layer.  It  is  shown  in  this  figure  that  the  velocity  of 
the  flow  near  the  wall  decreased  rapidly  with  Increasing  distance  along  the  model 
for  values  of  local  surface  angle,  as ,  between  0  and  l8  deg.  However,  the  velocities 
near  the  wall  decreased  very  little  between  values  of  cs  of  l8  and  30  deg  because  of 
increased  mixing  of  the  flow  into  the  inner  portion  of  the  boundary  layer. 

The  effect  of  changes  in  the  boundary  layer  thickness  at  Station  0  on  boundary 
layer  profiles  at  the  downstream  end  of  the  curved- surf ace  portion  of  the  models  is 
given  in  Fig  23.  An  increase  In  boundary  layer  thickness  at  Station  0  resulted  in 
an  Increase  in  boundary  layer  thickness  and  a  distortion  of  profile  shape  at  the 
downstream  end  of  the  model.  Two  of  the  profiles  shown  in  Fig.  23  were  obtained 
during  tests  in  which  the  boundary  layer  was  separated  from  the  surface  of  the  model 
near  the  station  where  the  local  model  surface  angle  was  15  deg.  These  two  profiles 
were  not  included  in  the  series  of  38  zero-bleed  runs  listed  in  Table  I. 

The  variations  of  boundary  layer  displacement  thickness,  8*,  boundary  layer 
momentum  thickness,  9  ,  and  reduced  boundary  layer  shape  parameter,  G  ,  with  distance 
from  the  curved- surface  tangent  point  are  shown  in  Fig.  2k  for  three  different  values 
of  initial  boundary  layer  thickness.  The  reduced  boundary  layer  shape  parameter,  G, 
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increased  with  increasing  distance  along  the  wall  for  values  of  local  surface  angle, 
as,  between  approximately  0  and  15  deg  (approximately  the  same  range  of  as  for  which 
the  velocities  near  the  wall  decreased  in  Fig.  22)  and  decreased  for  values  of  as 
greater  than  approximately  15  deg  (where  the  velocities  near  the  wall  in  Fig.  22 
were  relatively  uninfluenced  by  distance  along  the  curved  surface).  The  significance 
of  this  observation  is  discussed  in  a  following  section. 

Typical  schlieren  photographs  obtained  during  tests  of  the  two-radii  two- 
dimensional  curved- surface  models  are  shown  in  Fig.  25-  The  variations  of  Mach  number, 
boundary  layer  momentum  thickness,  9  ,  and  reduced  boundary  layer  shape  parameter,  G, 
with  distance  along  the  wall  of  the  two-radii  models  are  given  in  Fig.  26. 

A  schlieren  photograph  obtained  during  tests  of  the  axisymmetric  curved-surface 
models  is  shown  in  Fig.  27.  The  wave  patterns  are  similar  to  those  shown  in 
Figs.  19  and  25  for  the  two-dimensional  models.  The  variations  of  wall  static  pres¬ 
sure  and  Mach  number  with  distance  from  the  curved- surface  tangent  point  of  a 
typical  axisymmetric  curved- surface  model  are  given  in  Fig.  28.  The  theoretical 
variations  of  static  pressure  and  Mach  number  with  distance  shown  therein  were  ob¬ 
tained  from  a  method-of- characteristics  solution  of  the  axisymmetric  inviscid  flow 
over  the  model.  The  variations  of  boundary  layer  displacement  thickness,  8*, 
boundary  layer  momentum  thickness,  6  ,  and  reduced  boundary  layer  shape  parameter,  G, 
with  distance  from  the  curved- surface  tangent  point  of  an  axisymmetric  curved- 
surface  model  are  shown  in  Fig.  29.  The  increase  in  local  radius  (in  a  plane  normal 
to  the  axis  of  symmetry)  with  distance  along  the  axisymmetric  curved- surface  models 
should  result  in  a  thinning  out  of  the  boundary  layer  on  the  curved  surface  if  the 
boundary  layer  displacement  "area"  and  momentum  "area"  are  assumed  to  remain  constant. 
(As  noted  in  the  section  entitled  ANALYTICAL  INVESTIGATION,  it  is  possible  for 
boundary  layer  momentum  and  displacement  thickness  areas  to  decrease  in  an  adverse 
pressure  gradient.)  The  dashed  lines  in  Fig.  29,  which  were  obtained  by  multiplying 
the  initial  values  of  momentum  thickness  and  displacement  thickness  by  the  ratio  of 
initial  radius  to  local  radius,  indicate  whether  the  momentum  area  and  displacement 
area  increased  or  decreased  along  the  length  of  the  model. 


Boundary  Layer  Characteristics  With  Bleed 

A  series  of  tests  was  conducted  to  determine  the  effect  of  boundary  layer  removal 
on  the  boundary  layer  characteristics  of  some  of  the  single-radius  two-dimensional 
curved- surface  models.  The  models  employed  had  radii  of  1.5  and  3  in.  (see  Figs.  3 
and  b) .  Tests  were  conducted  with  different  values  of  approach  boundary  layer 
thickness  (as  determined  by  the  pressure  in  the  injection  plenum)  and  with  different 
amounts  of  bleed  (as  determined  by  the  number  and  location  of  open  bleed  holes  in 
the  model).  The  number  of  holes  employed  in  a  given  test  was  controlled  by  filling 
a  fraction  of  the  total  number  of  holes  with  a  plastic  compound. 
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Typical  schlieren  photographs  of  the  boundary  layer  bleed  models  are  shown  in 
Fig.  30.  These  photographs  show  the  flow  pattern  with  no  boundary  layer  bleed 
(Fig.  30a),  the  flow  pattern  with  insufficient  bleed  to  prevent  separation  (Fig.  30b), 
and  the  flow  pattern  with  just  enough  boundary  layer  bleed  to  prevent  separation 
(Fig.  30c).  The  exact  combination  of  values  of  initial  boundary  layer  thickness  and 
bleed  flow  which  caused  incipient  separation  was  determined  by  observation  of  the 
behavior  of  an  oil  film  which  was  injected  through  one  of  the  wall  static  pressure 
taps.  Incipient  separation  was  indicated  by  an  accumulation  of  oil  on  the  model 
surface  at  the  separation  point. 

The  schlieren  photograph  in  Fig.  31  shows  the  flow  pattern  when  the  bleed  flow  was 
considerably  greater  than  that  necessary  to  prevent  separation.  Under  these  conditions 
a  shock  wave  originated  from  each  set  of  bleed  holes.  No  waves  were  noted  in  any  tests 
in  which  the  bleed  flow  was  equal  to  or  less  than  that  for  incipient  separation. 

Preliminary  tests  indicated  that  the  amount  of  bleed  flow  required  to  prevent 
separation  was  less  if  the  bleed  holes  were  located  between  local  surface  angles  of 
12  and  24  deg  than  if  the  bleed  holes  were  located  at  local  surface  angles  less  than 
12  deg  or  more  than  2b  deg.  Therefore,  all  data  with  boundary  layer  bleed  presented 
in  this  report  were  obtained  with  the  bleed  holes  located  between  local  surface 
angles  of  12  and  2b  deg. 

Typical  wall  static  pressure  and  Mach  number  distributions  along  a  two- 
dimensional  curved- surface  model  employing  boundary  layer  bleed  are  given  in  Fig.  32- 
Typical  boundary  layer  characteristics  are  shown  in  Figs.  33  and  34.  The  data  pre¬ 
sented  in  Figs.  32  througn  3^-  were  obtained  with  a  bleed  rate  approximately  twice 
that  necessary  to  prevent  separation  on  this  model.  It  was  not  possible  to  obtain 
boundary  layer  data  at  all  stations  on  the  model  at  the  incipient-separation  condi¬ 
tion  because  the  presence  of  the  probe  caused  local  separation  of  the  boundary  layer. 

The  effect  on  boundary  layer  velocity  profiles  of  bleeding  more  boundary  layer 
from  the  model  than  the  amount  required  to  prevent  separation  is  shown  in  Fig.  35  for 
a  local  surface  angle  of  24  deg.  This  figure  shows  that  the  use  of  more  boundary 
layer  bleed  than  that  required  to  prevent  separation  had  very  little  effect  on  the 
values  of  reduced  boundary  layer  shape  parameter,  G  ,  at  this  station. 


ANALYTICAL  INVESTIGATION 


A  complete  determination  of  the  characteristics  of  a  boundary  layer  in  an  adverse 
pressure  gradient  requires  knowledge  of  both  the  absolute  thickness  of  the  boundary 
layer  (as  measured  by  the  boundary  layer  momentum  thickness,  6  )  and  of  the  shape  of 
the  boundary  layer  (as  measured  by  the  boundary  layer  shape  parameter,  H  ,  or  the  re¬ 
duced  boundary  layer  shape  parameter,  G  -  see  Eq.  l).  Although  the  problems  of  de¬ 
termining  8  and  G  are  not  independent  of  each  other,  they  are  treated  separately  in 
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the  following  two  subsections.  It  is  assumed  in  the  discussion  on  the  following  pages 
that  the  total  temperature  in  the  boundary  layer  is  constant.  A  discussion  of  the 
effect  of  wall  temperature  on  boundary  layer  characteristics  is  given  in  Appendix  I. 


Calculation  of  Boundary  Layer  Momentum  Thickness 

A  conventional  form  of  the  von  Karman  boundary  layer  momentum  equation  is  given 
in  Ref.  31  as  follows: 


cfq  =  -^(pvze)  +  pvs*-^. 


(2) 


Rearrangement  of  this  equation  leads  to  the  following: 


d0 


cf  ,  _  rdq/q  dv/v 

d*_erdsr  +  H“dir 
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(3) 


The  first  term  on  the  right  side  of  the  equation  represents  the  change  in  8  due  to  a 
change  in  distance,  whereas  the  second  term  on  the  right  side  of  the  equation  repre¬ 
sents  the  change  in  8  due  to  a  change  in  Mach  number.  The  factors  affecting  the 
two  right-side  terms  will  be  treated  separately  in  the  following  discussion. 

The  local  skin-friction  coefficient  in  a  region  of  zero  Mach  number  gradient  is 
given  in  Fig.  3 6  as  a  function  of  Mach  number  and  the  Reynolds  number  based  on  the 
momentum  thickness  of  the  boundary  layer.  The  variations  of  local  skin-friction 
coefficient  with  Mach  number  and  with  Reynolds  number  shown  in  this  figure  were  ob¬ 
tained  by  cross  plotting  the  experimental  data  from  Refs.  33  through  36.  The  curves 
of  friction  coefficient  shown  for  the  transition  region  between  laminar  flow  and 
turbulent  flow  were  drawn  arbitrarily  through  the  small  number  of  data  points  which 
were  available . 

For  incompressible  flow,  it  is  known  that  the  friction  coefficient  in  a  region 
of  adverse  pressure  gradient  is  considerably  less  than  the  friction  coefficient  on 
a  flat  plate  at  the  same  Reynolds  number  because  of  the  decrease  in  velocity  gradient 
at  the  wall  caused  by  the  adverse  pressure  gradient.  According  to  Ref.  37  j  the 
variation  of  friction  coefficient  with  H  for  incompressible  flow  is  as  follows: 
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cf 


I 

6'.54h 


(4) 


Rearranging  yields 


cf 


i  .  i  . _ !_ 

el.56HN=7G  "  el.56(l.287)G  ‘  g2.0lG 


(5) 


To  account  for  the  fact  that  G  on  a  flat  plate  is  usually  different  from  unity  (see 
Fig.  l8),  Eq.  (5)  has  been  changed  arbitrarily  to  the  following  form: 


cf  w  2.01675? 


(6) 


This  equation  is  plotted  in  Fig.  37-  The  curves  in  Fig.  37  for  Mach  numbers  other 
than  zero  were  determined  in  a  comparable  manner  by  the  use  of  the  transformation 
employed  in  Ref.  26. 

As  noted  in  a  preceding  paragraph,  the  second  term  on  the  right  side  of  Eq.  (3) 
represents  the  change  in  9  due  to  a  change  in  Mach  number.  This  term  can  be  sim¬ 
plified  by  the  introduction  of  a  quantity  called  the  boundary  layer  deceleration 
parameter,  J ,  where 


J  = 


dq/q  ,  ..dv/v 

dM  dM 


(7) 


It  can  be  shown  that,  for  isentropic  flow  outside  the  boundary  layer. 


dq/q 

dM 


2_ 

M 


r  M 


1  + 


r- 1 


M 


dwv..  j_  _  r 2  M 

dM  -  M  |+2L|L  m2 


(8) 


(9) 
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Substituting  Eqs.  (l),  (8),  and  (9)  into  Eq.  (7)  yields 


For  a  value  of  the  ratio  of  specific  heats,  y  =  1.4,  Eq.  (10)  becomes 


The  boundary  layer  deceleration  parameter,  J,  as  evaluated  herein,  is  given  in 
Figs.  38  through  4l  as  a  function  of  Mach  number  and  the  reduced  boundary  layer  shape 
parameter,  G. 

Substitution  of  Eq.  (7)  in  Eq.  (3)  and  conversion  from  differential  to  incre¬ 
mental  form  yields 

A9=  ^-AX-  6 JAM  (12) 

This  equation  can  b  used  in  a  step-wise  calculation  of  the  boundary  layer  growth  in 
an  adverse  Mach  number  gradient.  However,  in  evaluating  this  relation,  if  the 
momentum  thickness,  9  ,  at  the  beginning  of  each  step  is  used  rather  than  an  average 
value  of  9  during  the  step,  errors  are  introduced  into  the  calculation  These  errors 
can  be  avoided  by  the  use  of  an  additional  factor  which  is  developed  in  Appendix  II. 
This  additional  factor,  which  is  a  correction  factor  for  the  variation  of  9  during 
the  step,  converts  Eq.  (12)  to  the  following  form: 


A0*[^-  AX-0,JAm][|-~^]  (13) 
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The  variation  of  boundary  layer  momentum  thickness  with  Mach  number  for  zero 
wall  friction  is  given  in  Figs.  42,  4-3,  and  44  for  three  different  values  of  G 
(assumed  constant  during  the  deceleration  processes)  and  for  three  different  values 
of  initial  Mach  number.  The  data  shown  in  these  figures  were  obtained  by  a  step-wise 
integration  using  Eq .  (13).  The  results  for  a  value  of  G  of  unity  were  checked  by 
the  use  of  the  tables  in  Ref.  31  neglecting  the  friction  term  in  that  reference. 

These  three  figures  show  that,  for  zero  wall  friction  and  a  value  of  G  of  unity,  a 
reduction  in  Mach  number  results  in  a  decrease  in  momentum  thickness  for  Mach  numbers 
above  2.5  and  an  increase  in  momentum  thickness  for  Mach  numbers  below  2.5.  This 
same  result  could  be  obtained  from  Fig.  38,  where  the  boundary  layer  deceleration 
parameter,  J,  changes  sign  at  a  Mach  number  of  2.5  for  a  value  of  G  of  unity. 

Figure  38  shows  that  an  increase  in  G  results  in  an  increase  in  the  Mach  number  at 
which  J  equals  zero,  i.e.,  an  increase  in  the  Mach  number  for  which  a  decrease  in 
Mach  number  results  in  no  change  in  boundary  layer  momentum  thickness. 


Calculation  of  Boundary  Layer  Profile  Shape 
Conventional  Procedures 


The  boundary  layer  profile  shape  parameter,  H  (or  the  reduced  shape  parameter, 
G)  must  be  known  for  three  reasons: 

1)  H  influences  the  growth  of  Q  in  an  adverse  pressure  gradient  (see  preceding 
subsection) . 

2)  H,  in  combination  with  Q,  determines  the  boundary  layer  displacement  thick¬ 
ness,  8*,  which  is  needed  to  evaluate  the  amount  that  a  wall  must  be 
diverged  to  account  for  boundary  layer  growth  so  that  the  inviscid  flow 
outside  the  boundary  layer  is  unchanged  by  the  boundary  layer  growth. 

3)  H,  relative  to  H  for  a  flat  plate  (see  Figs.  13  and  18),  determines  whether 
or  not  the  boundary  layer  is  close  to  being  separated  from  the  wall. 

The  boundary  layer  shape  parameters,  H  and  G,  are  functions  of  Reynolds  number  and 
the  history  of  the  wall  pressure  gradient.  The  effect  of  Reynolds  number  on  the 
value  of  G  for  a  flat  plate  with  zero  pressure  gradient  is  given  in  Fig.  l8. 

The  effect  of  adverse  pressure  gradients  on  the  boundary  layer  shape  parameter, 
H,  for  incompressible  flow  is  usually  calculated  from  an  equation  having  the  follow¬ 
ing  form: 


dH*f,  (H)^-  +f2(H,Rg)^- 


(14) 
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The  first  term  on  the  right  side  of  Eq.  ( l4 )  results  in  an  increase  in  H  for  a  decrease 
in  velocity.  The  second  term  on  the  right  side  results  in  a  decrease  in  H  for  an  in¬ 
crease  in  X  (i.e.,  the  second  term  is  usually  negative).  Different  forms  of  the  first 
and  second  terms  on  the  right  side  are  recommended  by  different  investigators  (see 
Refs.  27  through  30).  The  particular  forms  of  Eq.  (l4)  used  in  Refs.  27  and  30, 
respectively,  for  the  insulated-wall  case  follow: 

dH  =  -  HiH+l)2  (H~  ^  -0.00369  Rl0268  — 

2  v  e  56  H  0  9 


dH*  -e5(H_l4)  Q  -o.oi35e5(H''-4)  (h-i.4)  rJ 


(15) 


(16) 


The  form  of  Eq.  (14)  employed  in  Ref.  27,  Eq.  (15),  is  extended  in  Ref.  26  to  com¬ 
pressible  flow  by  a  transformation  procedure. 

Lag-Length  Procedure 

The  lag-length  procedure  discussed  in  th.  following  paragraphs  for  calculating 
G  in  adverse  pressure  gradients  was  developed  because  of  difficulties  encountered  in 
the  use  of  Eq.  (l4)  -  see  the  following  section.  Development  of  this  procedure  re¬ 
quired  knowledge  of  the  amount  of  flow  in  the  boundary  layer  in  the  adverse  pressure 
gradient.  The  variation  with  Mach  number  of  the  total  flow  within  a  boundary  layer 
for  initial  Mach  numbers  of  3*0,  5-0,  and  10.0  is  given  in  Fig.  45  for  a  constant 
value  of  G  of  unity  during  the  deceleration  process.  The  amount  of  flow  within  the 
boundary  layer  was  evaluated  from  the  values  of  9  given  in  Figs.  42  through  44  and 
from  the  tables  of  boundary  layer  parameters  for  N  =  7  in  Ref.  31-  The  increase  of 
flow  in  the  boundary  layer  with  a  decrease  in  Mach  number  shown  in  this  figure  repre¬ 
sents  the  influx  of  flow  necessary  to  increase  the  average  total  pressure  in  the 
boundary  layer  from  the  low  value  associated  with  a  high  Mach  number  G  =  1  profile 
to  the  high  value  associated  with  a  low  Mach  number  G  =  1  profile.  Since  the 
variation  of  boundary  layer  momentum  thickness,  9  ,  with  Mach  number  in  Figs.  42 
through  44  was  calculated  on  the  basis  of  zero  wall  friction,  the  variation  of  total 
flow  in  the  boundary  layer  with  Mach  number  given  in  Fig.  45  does  not  include  the 
influx  of  flow  into  the  boundary  layer  necessary  to  overcome  the  shear  force  at  the 
wall . 


An  example  of  the  actual  variation  of  flow  within  the  boundary  layer  with  dis¬ 
tance  along  the  surface  of  a  curved- surface  model  is  given  in  Fig.  46.  Two  reference 
curves  as  well  as  the  measured  curve  are  given  in  this  figure:  one  shows  the  varia- 
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tion  of  flow  within  the  boundary  layer  with  distance  which  would  be  required  to 
maintain  a  G  =  Gf  profile  shape  with  zero  pressure  gradient;  the  second  reference 
curve  shows  the  variation  with  distance  of  the  flow  within  the  boundary  layer  which 
would  be  required  to  maintain  a  G  =  Gp  profile  for  the  measured  pressure  gradient. 

This  figure  shows  that  the  variation  of  the  actual  flow  within  the  boundary  layer 
with  distance  along  the  curved  surface  is  approximately  equal  to  that  expected  on  a 
flat  plate  for  values  of  (X-Xo)  between  0  and  1-5  in.,  and  that  the  amount  of  flow 
within  the  boundary  layer  increases  rapidly  only  for  values  of  (X-X0)  larger  than 
1.5  in.  It  can  be  seen  from  Figs.  24  and  46  that  the  region  where  the  rate  of  influx 
of  flow  into  the  boundary  layer  in  the  latter  figure  is  small,  (X-Xo)  from  0  to 
1.5  in.,  corresponds  to  a  region  of  rapidly  increasing  G  in  the  former  figure  and 
that  the  region  in  which  the  influx  of  flow  into  the  boundary  layer  in  Fig.  46  is 
large  (downstream  of  (X-Xq)  =  1.5)  corresponds  to  a  region  of  constant  or  decreasing 
G  in  Fig.  24. 

It  is  possible  to  explain  why  the  position  at  which  flow  actually  starts  to 
enter  the  boundary  layer  lags  behind  the  position  at  which  flow  should  come  into  the 
boundary  layer  to  maintain  a  G  =  Gp  profile  as  follows:  the  adverse  pressure 
gradient  is  first  felt  in  the  boundary  layer  as  a  decrease  in  velocity  near  the  wall 
which  is  greater  (percentage -wise)  than  the  decrease  in  velocity  of  the  flow  near 
the  outer  edge  of  the  boundary  layer  (see  Fig.  22).  It  Is  hypothesized  that  the  re¬ 
sulting  change  in  velocity  near  the  wall  from  the  velocity  that  would  exist  In  a 
flat-plate  velocity  profile  induces  a  turbulence  in  order  to  cause  the  high-energy 
flow  near  the  outer  part  of  the  boundary  layer  to  mix  with  the  low-energy  flow  near 
the  wall.  It  is  also  hypothesized  that  some  distance  Is  required  for  this  increased 
turbulence  to  move  to  the  outer  edge  of  the  boundary  layer  and  that  this  distance  is 
equal  to  1.5  in.  for  the  data  of  Fig.  46.  Some  data  to  support  this  hypothesis  are 
given  in  Ref.  38,  where  it  is  shown  that  the  height  of  the  region  of  high  turbulence 
in  the  boundary  layer  of  a  subsonic  diffuser  increases  relative  to  the  boundary 
layer  thickness  as  the  distance  from  the  entrance  of  the  diffuser  is  increased. 

A  tentative  assumption,  suggested  in  Ref.  39>  is  that  the  rate  of  mixing  of  flow 
into  the  boundary  layer  in  an  adverse  pressure  gradient  is  approximately  equal  to 
that  on  a  flat  plate.  It  appears  from  Fig.  46  that  this  assumption  is  accurate  only 
near  the  upstream  end  of  the  adverse  pressure  gradient  (corresponding  to  (X-Xq)  be¬ 
tween  0  and  1.5  in.)  and  that  this  assumption  is  extremely  inaccurate  near  the  down¬ 
stream  end  of  the  adverse  pressure  gradient. 

Additional  plots  of  the  flow  within  the  boundary  layer,  prepared  from  the  experi¬ 
mental  data  presented  in  this  report,  showed  the  same  trend  as  that  shown  in  Fig.  46 
but,  in  some  cases,  exhibited  considerable  scatter  because  of  inaccuracies  associated 
with  the  definition  of  the  exact  edge  of  the  boundary  layer.  Consequently,  an 
alternate  method  of  consistently  defining  the  flow  within  the  boundary  layer  was 
needed  which  did  not  require  definition  of  the  edge  of  the  boundary  layer.  It  is 


18 


R-1285-U 


possible  to  obtain  a  consistent  measure  of  the  flow  within  a  boundary  layer  from  8* 
and  6  if  it  is  assumed  that  all  boundary  layer  profiles  are  members  of  the  same 
family  (it  is  conventional  in  theoretical  analyses  of  boundary  layers  to  make  the 
assumption  of  a  common  profile  family).  The  profile  family  adopted  in  this  report 
is  assumed  to  be  made  up  of  profiles  similar  to  the  equivalent  step  profile  shown 
in  Fig.  kj ,  As  noted  therein,  the  equivalent  step  profile  has  the  same  values  of 
8*  and  8  and,  hence,  the  same  values  of  H  and  G  as  the  actual  boundary  layer  profile. 
However  the  total  boundary  layer  thickness  and  the  total  flow  within  the  equivalent 
step  profile  are  less  than  those  for  the  actual  boundary  layer  profile.  Values  of 
the  ratio  of  the  velocity  within  equivalent  step  profiles  to  the  velocity  outside 
the  bounaary  layer  are  given  in  Fig  48  for  different  values  of  Mach  number  outside 
the  boundary  layer.  The  equivalent  step  profile  was  employed  as  a  tool  in  this 
report  because  it  is  easily  adaptable  to  physical  reasoning  concerning  the  behavior 
of  boundary  layers  in  adverse  pressure  gradients  A  discussion  of  the  influence  of 
employing  a  family  of  profiles  other  than  that  for  equivalent  step  profiles  is  given 
in  the  following  paragraphs . 

It  is  convenient  in  the  development  of  the  lag-length  procedure  to  first  con¬ 
sider  the  behavior  of  an  equivalent  step  profile  in  a  Mach  number  gradient  corres¬ 
ponding  to  the  region  of  (X-Xo)  between  0  and  1.5  in.  in  Fig.  46,  i.e.,  in  a  region 
in  which  the  deceleration  is  so  rapid  that  the  rate  of  mixing  of  flow  into  the 
boundary  layer  from  the  free  stream  is  not  increased  above  that  which  occurs  upstream 
of  the  adverse  pressure  gradient-  Since  the  influx  of  flow  into  the  actual  boundary 
layer  in  this  region  is  approximately  equal  to  that  necessary  to  overcome  the  shear 
force  at  the  wail,  it  is  permissible  to  assume  that  there  is  no  influx  of  flow  into 
the  equivalent  step  profile  and  no  shear  force  at  the  wall  for  the  equivalent  step 
profile  in  this  region.  This  assumption  is  equivalent  to  stating  that  the  ratio  of 
the  total  pressure  within  the  equivalent  step  profile  to  the  total  pressure  outside 
of  the  equivalent  step  profile  is  constant  during  the  deceleration  process.  However, 
it  can  be  shown  from  isentropic  flow  relations  that  the  ratio  of  the  velocity  within 
the  equivalent  step  profile  to  the  velocity  outside  the  boundary  layer  decreases 
rapidly  during  such  a  deceleration  process.  (The  velocity  in  the  equivalent  step 
profile  approaches  zero  when  the  flow  has  decelerated  to  the  condition  where  the 
wall  static  pressure  is  equai  to  the  total  pressure  within  the  equivalent  step 
profile .  ) 

Tne  variations  of  the  boundary  layer  shape  parameters,  H  and  G,  with  fractional 
decrease  in  velocity  for  incompressible  fxow  calculated  on  the  basis  of  a  constant 
ratio  of  step -prof ile  total  pressure  to  total  pressure  outside  the  boundary  layer 
are  given  by  the  so. id  line  in  Fig  49-  Also  shown  in  this  figure  are  the  variations 
of  H  and  G  with  fractional  decrease  in  velocity  calculated  from  the  first  term  on 
the  right  side  of  Eqe .  (15)  and  (l6)  The  figure  shows  that  the  variations  of  H  and 
G  with  fractional  decrease  in  Mach  number  are  of  the  same  general  form  for  the 
equivaxent  step  profixe  as  those  calculated  from  the  first  term  on  the  right  side  of 
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Eqs.  (15)  and.  (16).  Therefore,  the  first  term  of  Eqs.  (1*0,  (15),  or  (16),  which  is 
much  larger  than  the  second  terms  of  these  equations  for  rapid  changes  in  Mach  number 
(small  AX),  represents  the  change  in  boundary  layer  shape  which  occurs  when  the  rate 
of  influx  of  flow  into  the  boundary  layer  has  not  been  increased  above  that  which 
would  be  obtained  on  a  flat  plate.  The  effect  of  a  rapid  decrease  in  Mach  number  on 
the  reduced  boundary  layer  shape  parameter,  G,  for  compressible  flow  is  shown  in 
Fig.  50.  The  variation  of  G  with  fractional  decrease  in  Mach  number  shown  in  Figs.  49 
and  50  applies  only  when  the  decrease  in  Mach  number  is  so  rapid  that  the  mixing  which 
occurs  downstream  of  (X-Xq)  of  1.5  in.  in  Fig.  46  has  no  influence  on  the  profile. 

In  order  to  use  the  characteristics  of  equivalent  step  profiles  in  regions  where 
there  is  a  large  influx  of  flow  into  the  boundary  layer  (such  as  the  region  down¬ 
stream  of  (X-Xq)  of  1.5  in.  in  Fig.  46),  a  new  boundary  layer  parameter  was  adopted. 
This  new  parameter,  called  the  boundary  layer  flow  thickness,  (3 ,  is  related  to  the 
flow  within  an  equivalent  step  profile  by  the  following  equation: 


0.53IPt  B 

Wesp'“/?7^ 


(17) 


The  total  pressure  and  total  temperature  in  the  preceding  equation  are  the  total 
pressure  and  total  temperature  outside  the  boundary  layer.  Therefore,  (3  represents 
the  height  of  a  sonic  nozzle  which  will  pass  the  flow  within  an  equivalent  step 
profile  at  free-stream  total  pressure  and  total  temperature.  The  variations  of  di¬ 
mensionless  flow  thickness,  (3/8,  with  G  for  various  Mach  numbers  outside  the  boundary 
layer  are  given  in  Figs.  51  and  52.  The  curves  on  these  figures  were  calculated  for 
a  value  of  X  of  1.4  and  a  value  of  total  temperature  within  the  step  profile  equal 
to  the  total  temperature  outside  the  step  profile.  The  flow  within  an  equivalent 
step  profile  (as  measured  by  (3)  at  a  given  Mach  number  can  be  determined  from  the 
following  equation: 


(18) 


An  example  of  the  variation  of  (3  with  distance  from  the  curved- surface  tangent 
point  is  given  in  Fig.  53*  (The  data  employed  to  construct  this  figure  are  the  same 
as  those  employed  for  Fig.  46.)  The  method  of  determination  of  (3  in  Fig.  53  is  noted 
on  the  figure.  It  can  be  seen  that  the  relationship  between  the  actual  and  ideal 
values  of  /3  shown  therein  is  qualitatively  the  same  as  that  for  the  corresponding 
curves  in  Fig.  46.  The  displacement  of  tne  actual  /3  curve  from  the  ideal  (3  curve. 
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shown  on  Fig.  53,  is  described  in  this  report  in  terms  of  the  lag  length, The 
method  of  determination  of  ^  from  the  experimental  data  is  noted  on  the  figure.  The 
lines  between  the  actual  0  and  ideal  0  curves  on  Fig.  53  are  drawn  at  a  slope  which 
is  calculated  from  the  rate  of  change  of  ^  on  a  flat  plate  at  the  beginning  of  the 
adverse  pressure  gradient.  The  values  of  lag  length, ,  determined  from  this  con¬ 
struction  are  shown  in  the  lower  part  of  the  figure  in  a  nondimensionalized  form. 

The  lag  lengths  are  nondimensionalized  by  dividing  by  the  momentum  thickness  at  the 
upstream  end  of  the  lag  length. 

The  experimental  values  of  G  shown  in  Fig.  29  can  be  back- calculated  from  the 
lag-length  data  shown  in  Fig.  53  by  the  inverse  of  the  procedure  by  which  the  lag 
lengths  were  obtained.  Therefore,  the  generalized  problem  of  determining  G  in  an 
adverse  pressure  gradient  is  reduced  to  one  of  determining  the  dimensionless  lag 
length,  Ji/B  . 

A  simplified  procedure  for  approximating  the  boundary  layer  shape  parameter,  G, 
in  small  adverse  pressure  gradients  can  also  be  derived  on  the  basis  of  lag-length 
data.  This  procedure,  which  is  illustrated  in  Fig.  54,  is  based  on  the  use  of  curves 
such  as  those  presented  in  Figs.  49  and  50.  The  value  of  G  at  a  point  such  as  "b"  in 
Fig.  54  is  approximately  equal  to  the  value  of  G  which  would  be  obtained  by  rapid 
deceleration  of  the  flow  from  a  flat-plate  profile  ~.t  Point  "a".  The  difference  in 
Mach  numbers  between  Points  "a"  and  "b"  can  then  be  used  in  curves  such  as  those 
presented  in  Figs.  49  or  50.  The  difference  in  axial  distance  between  Points  "a" 
and  "b"  is  given  by  the  lag  length. 

Although  the  equivalent  step  profile  shown  in  Fig.  47  was  used  extensively  in 
the  development  of  the  reasoning  behind  the  lag-length  procedure,  the  only  place 
where  the  characteristics  of  the  equivalent  step  profile  are  actually  used  in  de¬ 
termining  lag  length  is  through  the  curves  of  0/8 given  in  Figs.  51  and  52.  It  is 
possible  to  generate  curves  similar  to  those  shown  in  Figs  51  and  52  for  a  different 
family  of  profile  shapes  and,  hence,  to  generate  lag-length  data  which  are  different 
from  those  shown  in  Fig.  53 •  However,  if  a  new  family  of  profile  shapes  were  selected 
in  which 0/8  was  a  constant  times  the  values  shown  in  Figs.  51  and  52,  then  the  lag 
lengths  determined  by  a  procedure  similar  to  that  described  in  Fig.  53  would  be 
identical  to  those  determined  from  the  equivalent  step  profile  because  of  the  method 
employed  to  attain  lag  lengths  from  the  0/8  data  (see  Fig.  53).  Evaluation  of  the 
flow  within  power  profiles  from  the  tables  of  Ref.  31  indicates  that  the  flow  within 
a  power  profile  for  a  given  value  of  6  (a  measure  of  /3/d)  is  approximately  twice  that 
for  an  equivalent  step  profile  for  a  wide  range  of  Mach  numbers  and  for  values  of 
N  between  5  and  11.  Lag  lengths  different  from  those  determined  from  an  equivalent 
step  profile  would  result  only  from  the  use  of  different  types  of  variations  of  0 /B 
with  Mach  number  and  G  from  those  shown  in  Figs.  51  and  52.  Although  values  of  JL/8 
for  the  data  presented  in  this  report  have  not  been  evaluated  for  any  variation  of 
0/8  with  Mach  number  and  G  other  than  that  given  in  Figs.  51  and  52  because  of  the 
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large  amount  of  labor  Involved,  a  cursory  examination  of  the  use  of  alternate  profile 
shapes  for  the  evaluation  of /3/S  has  Indicated  that  the  scatter  of  the  data  on  ■£■/ 9 
presented  in  this  report  would  not  be  reduced  significantly  by  the  use  of  alternate 
profile  shapes  for  the  evaluation  of  /3/S.  It  is  believed  that  significant  improve¬ 
ments  in  the  correlation  of  data  from  the  correlation  based  on  /S/S  presented  in  this 
report  can  only  be  made  by  the  use  of  multi-family  profile  shapes  (such  as  a  multi- 
step  profile  with  no  restriction  on  the  relative  values  of  the  velocities  in  each 
step).  A  correlation  of  data  based  on  the  use  of  multi-family  profile  shapes  would 
be  extremely  complicated. 


CORRELATION  OF  EXPERIMENTAL  AND  ANALYTICAL  RESULTS 


Boundary  Layer  Profile  Shape 

Determination  of  boundary  layer  profile  shape  for  an  arbitrary  configuration  by 
the  lag-length  procedure  discussed  in  the  preceding  section  requires  knowledge  of 
the  lag-length  characteristics  of  the  boundary  layer  for  the  arbitrary  configuration. 
Data  on  lag  length  have  been  obtained  from  the  38  runs  without  boundary  layer  bleed 
(see  Table  I)  in  order  to  obtain  a  correlation  of  lag  length  which  can  be  used  for 
an  arbitrary  configuration.  Typical  values  of  boundary  layer  flow  thickness  from 
which  lag  length  can  be  determined  (see  Fig.  53)  and  lag  length  are  presented  in 
Figs.  55  and  56  for  a  single-radius  two-dimensional  curved- surf ace  model,  in  Figs.  57 
and  58  for  a  two-radii  two-dimensional  curved- surface  model,  and  in  Fig.  59  for  an 
axisymmetric  curved- surface  model.  Although  the  scatter  in  the  values  of  lag  length 
for  a  given  run  was  large  as  a  result  of  the  scatter  in  the  values  of  G  which  were 
used  to  determine  lag  length  (see  Figs.  24,  26,  and  29),  it  appears  from  study  of 
the  data  from  all  38  runs  that  there  was  relatively  little  variation  of  lag  length 
with  distance  along  a  curved- surface  model  for  a  given  mn.  However,  there  appeared 
to  be  large  variations  of  lag  length  with  Reynolds  number  based  on  the  initial 
boundary  layer  momentum  thickness  (see  Figs.  56  and  58).  This  change  in  lag  length 
with  initial  Reynolds  number  was  particularly  evident  at  low  values  of  initial 
Reynolds  number .  It  is  believed  that  the  low  values  of  lag  length  at  low  values  of 
initial  Reynolds  number  are  caused  by  the  rapid  changes  in  boundary  layer  shape  in 
the  transition  region  which  occurs  at  Reynolds  numbers  based  on  momentum  thickness 
on  the  order  of  1000.  It  Is  hypothesized  that  the  vortices  which  cause  the  high 
mixing  rates  near  the  wall  to  propagate  to  the  outside  edge  of  the  boundary  layer 
(see  discussion  in  preceding  section)  are  much  larger  immediately  downstream  of 
transition  than  they  are  further  downstream  where  the  boundary  layer  has  a  fully 
developed  turbulent  profile.  These  large-size  vortices  would  require  less  lag  length 
than  small-size  vortices  to  mix  the  high-energy  flow  from  outside  the  boundary  layer 
with  the  relatively  low-energy  flow  near  the  wall. 
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Studies  of  various  methods  of  correlating  the  lag-length  data  have  indicated 
that  lag  length  is  primarily  a  function  of  the  Reynolds  number  at  the  beginning  of 
the  adverse  Mach  number  gradient  (see  Figs.  56  and  58)-  Although  no  data  were  ob¬ 
tained  in  which  there  was  an  extremely  large  variation  of  local  Reynolds  number 
during  a  test  with  a  single  model  at  a  single  Mach  number  and  initial  boundary  layer 
thickness,  it  was  hypothesized  that  the  influence  of  initial  Reynolds  number  should 
be  reduced  at  high  values  of  local  Reynolds  number.  This  hypothesis  must  be  true 
for  the  hypothetical  case  in  which  an  infinitesimally  small  pressure  gradient  is  im¬ 
posed  in  the  region  between  transition  (low  Rg)  and  the  beginning  of  a  severe  adverse 
pressure  gradient  at  high  Rg.  It  is  obvious  that  A/ 8  in  the  region  of  the  severe 
adverse  pressure  gradient  should  be  determined  from  evaluated  at  the  beginning  of 
the  severe  adverse  pressure  gradient  rather  than  at  transition.  Therefore,  in 
fitting  a  curve  through  the  lag-length  data,  it  was  decided  to  employ  an  equation  in 
which  the  effect  of  initial  Reynolds  numbers  would  die  out  at  high  values  of  local 
Reynolds  number.  The  equation  fitted  through  the  lag-length  data  obtained  from  the 
tests  described  in  this  report  is  as  follows: 


A  plot  of  the  preceding  equation  is  given  in  Fig.  60.  The  two  constants  in  the 
equation  (-22-5  and  +49-7)  were  determined  by  a  least-squares  fit  through  149  data 
points  obtained  from  the  38  runs  without  bleed  during  which  separation  did  not  occur. 
The  exponent  0.2  in  the  preceding  equation  was  arbitrarily  chosen  equal  to  the  ex¬ 
ponent  which  governs  the  variation  of  friction  coefficient  with  Reynolds  number  In 
turbulent  flow.  Calculations  in  which  the  exponent  0.2  was  replaced  by  smaller  and 
larger  values  indicated  that  the  accuracy  of  the  curve  fit  through  the  data  was  not 
changed  appreciably  for  values  of  this  exponent  between  0.1  and  0.4.  The  exponent 
0.5  in  the  preceding  equation  governs  the  rate  at  which  the  solid  curves  approach  the 
dashed  curve  in  Fig.  60  and  was  chosen  arbitrarily. 

A  comparison  of  lag  lengths  determined  from  the  data  and  from  the  curve  in  Fig.  6 0 
is  given  in  Fig.  6l  in  terms  of  the  lag-length  residual,^  .  As  noted  in  a  preceding 
paragraph,  the  scatter  of  the  lag-lengt’'  data  was  large  as  a  result  of  the  scatter  in 
the  values  of  G  employed  in  determining  lag  length  (see  Figs.  24,  26,  and  29).  The 
scatter  in  lag  length  as  indicated  by  the  lag-length  residuals  in  Fig.  6l  is  a  large 
percentage  of  the  absolute  value  of  lag  length  at  low  Reynolds  numbers  (see  Fig.  60). 
This  scatter  in  lag  length  at  low  Reynolds  numbers  results  from  the  large  change  in 
lag  length  with  Reynolds  number  at  low  Reynolds  number  and  the  difficulty  of  estab¬ 
lishing  the  exact  transition  point  upstream  of  an  adverse  pressure  gradient. 
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A  statistical  analysis  of  the  data  in  Fig.  6l  was  performed  to  determine  if  lag 
length  should  also  be  expressed  as  a  function  of  Mach  number,  type  of  model,  or 
dimensionless  pressure  gradient  (as  determined  by  the  ratio  of  average  model  radius 
to  initial  boundary  layer  momentum  thickness).  The  results  of  this  statistical 
analysis  are  given  in  Table  II.  The  average  residual  for  a  given  set  of  data 
in  the  table)  is  a  measure  of  how  much  the  average  of  this  set  of  data  falls  above 
or  below  the  fitted  curve  given  by  Eq.  (19) •  The  only  set  for  which  a  large  value 
of  average  residual  was  obtained  was  the  data  for  Mach  number  3-5>  The  large  average 
residual  for  the  Mach  number  3-5  data  is  not  believed  to  be  significant  because  of 
the  small  number  of  data  points  obtained  at  this  Mach  number  and  because  the  data  ob¬ 
tained  at  Mach  number  3-5  was  inconsistent  due  to  diff iculties  in  fixing  transition 
at  this  Mach  number.  Since  no  consistent  variation  of  average  lag-length  residual 
with  Mach  number,  model  type,  or  dimensionless  pressure  gradient  is  noted  in  the  data 
shown  in  Table  II,  it  is  concluded  that  none  of  these  variables  need  be  included  in 
the  general  expression  for  lag  length  given  by  Eq.  (19)* 

The  usefulness  of  the  lag-length  fairing  as  given  in  Eq.  (19)  and  in  Fig.  60  is 
determined  by  the  accuracy  with  which  the  relationship  for  lag  length  can  be  used  to 
calculate  the  variation  of  boundary  layer  shape  with  distance  in  an  adverse  pressure 
gradient.  Although  evaluation  of  G  by  the  lag-length  procedure  is  very  laborious  by 
hand  calculations,  such  an  evaluation  is  easily  adaptable  to  high-speed  computing 
machines.  Typical  comparisons  of  measured  values  of  reduced  boundary  layer  shape 
parameter  and  values  of  reduced  boundary  layer  shape  parameter  as  machine  calculated 
with  Eq.  (19)  are  given  in  Figs-  62,  63,  and  64  for  single-radius  two-dimensional, 
two-radii  two-dimensional,  and  axisymmetric  curved- surface  models,  respectively. 

Also  shown  in  these  figures  are  values  of  the  variation  of  reduced  boundary  layer 
shape  parameter  with  distance  along  the  curved  plates  as  calculated  by  the  method  of 
Ref.  26.  The  calculations  by  the  method  of  Ref.  26  produced  values cf  reduced  boundary 
layer  shape  parameter  greater  than  10  at  an  intermediate  point  in  the  adverse  pressure 
gradient  for  all  cases  given  in  Figs.  62  through  64.  The  method  of  Ref.  26,  which  is 
based  on  Eq.  (15),  gives  poor  answers  because  the  second  term  of  Eq.  (15)  does  not 
account  for  the  large  rates  of  mixing  of  flow  into  the  boundary  layer  near  the  down¬ 
stream  end  of  an  adverse  pressure  gradient.  The  agreement  between  the  measured 
values  of  G  and  the  values  of  G  calculated  from  the  lag-length  procedure  are  good, 
particularly  in  relation  to  the  large  differences  between  the  measured  values  of  G 
and  the  values  of  G  calculated  by  the  Ref.  26  method. 

A  summary  of  the  difference  between  the  maximum  measured  values  of  G  and  the 
maximum  values  of  G  calculated  by  the  lag-length  procedure  for  all  38  runs  for  which 
boundary  layer  bleed  was  not  employed  is  given  in  Fig.  65-  As  noted  thereon,  the 
average  difference  between  measured  and  calculated  values  of  maximum  G  is  0.044  and 
the  root-mean-square  difference  is  0.094.  If  the  data  for  Mach  number  3-5  were  omitted, 
the  average  difference  would  be  reduced  to  0.024  and  the  root-mean- square  difference 
would  be  reduced  to  O.O58. 
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The  calculated  values  of  reduced  boundary  layer  shape  parameter  shown  In  Figs.  62 
through  64  were  determined  on  the  basis  of  the  measured  variation  of  Mach  number  with 
distance  from  the  curved-surface  tangent  point  (see  Figs.  21,  26,  and  28).  The  effect 
of  small  changes  in  the  Mach  number  distribution  on  the  calculated  variation  of  the 
reduced  boundary  layer  shape  parameter  with  distance  is  shown  in  Fig.  66.  The  varia¬ 
tion  of  Mach  number  with  distance  given  by  Curve  1  of  this  figure  was  determined  from 
the  measured  variation  of  Mach  number  with  distance  given  in  Fig.  21.  The  variation 
of  Mach  number  with  distance  given  by  Curve  2  in  Fig.  66  was  determined  by  the 
theoretical  variation  of  Mach  number  with  distance  given  in  Fig.  21.  The  lower  plot 
in  Fig.  66  shows  that  the  small  difference  in  Mach  number  distribution  between  that 
given  by  Curve  1  and  that  given  by  Curve  2  results  in  an  appreciable  change  in  the 
variation  of  reduced  boundary  layer  shape  parameter  with  distance.  However,  the 
difference  between  the  maximum  value  of  reduced  boundary  layer  shape  parameter  for  the 
two  cases  is  only  0.05.  An  attempt  was  made  to  calculate  the  difference  between  the 
measured  and  theoretical  variations  of  Mach  number  with  distance  shown  in  Figs.  21 
and  66  on  the  basis  of  the  Mach  number  distribution  on  the  reference  flat  plate 
(Fig.  ll)  and  on  the  basis  of  the  Mach  number  changes  induced  by  the  variation  of 
boundary  layer  displacement  thickness  with  distance  along  the  curved  surface.  The 
resulting  corrected  theoretical  variation  of  Mach  number  with  distance  was  in  better 
agreement  with  the  measured  variation  of  Mach  number  with  distance  than  the  uncorrected 
theoretical  Mach  number  variation. 

Lag-length  data  were  also  evaluated  from  the  boundary  layer  profiles  given  in 
Refs.  9  and  l4.  The  boundary  layer  profiles  presented  in  Ref.  ik  were  obtained  for 
an  incident  oblique  shock  having  a  turning  angle  through  the  shock  of  10  deg  and  a 
Mach  number  of  2.97  upstream  of  the  shock.  The  Reynolds  number  based  on  the  boundary 
layer  momentum  thickness  approaching  the  shock  was  approximately  7400.  The  values  of 
-£/Q  determined  from  Ref.  l4  were  on  the  order  of  90,  which  is  almost  twice  the  value 
of  lag  length  for  the  correlation  given  in  Fig.  60.  The  difference  between  the  lag 
lengths  determined  from  Ref.  l4  and  those  of  this  report  might  be  due  to  the  fact 
that  the  Ref.  l4  variation  of  wall  static  pressure  with  distance  along  the  plate  was 
caused  by  shock  wave — boundary  layer  interaction,  whereas  the  present  variations  of 
wall  static  pressure  with  distance  along  the  curved  surface  were  caused  primarily  by 
the  local  slope  of  the  curved  surface. 

The  data  presented  in  Ref.  9  were  obtained  for  normal  shocks  with  a  Mach  number 
upstream  of  the  normal  shock  of  approximately  1-3-  The  Reynolds  numbers  based  on 
boundary  layer  momentum  thickness  approaching  the  shock  were  on  the  order  of  1200  to 
2300.  The  measured  values  ofX/Q  for  these  conditions  were  approximately  90,  as  com¬ 
pared  to  values  on  the  order  of  15  to  30  predicted  from  the  curve  presented  in  Fig.  60. 
No  satisfactory  explanation  for  this  difference  in  lag  length  has  been  devised  other 
than  that  suggested  for  the  Ref.  l4  data  discussed  in  the  preceding  paragraph. 
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Boundary  Layer  Momentum  Thickness 

According  to  the  boundary  layer  momentum  equation  (Eq.  3),  the  growth  of  boundary 
layer  momentum  thickness  in  an  adverse  pressure  gradient  is  a  function  of  the  wall 
friction  and  the  shape  of  the  boundary  layer  profile  during  the  deceleration  process. 
Five  sets  of  calculations  employing  different  assumptions  have  been  made  in  order  to 
evaluate  the  effect  of  changes  in  the  assumptions  regarding  wall  friction  and  boundary 
layer  profile  shape.  The  assumptions  employed  in  these  five  sets  of  calculations  are 
given  in  Fig.  67.  The  results  of  calculations  employing  the  five  setB  of  assumptions 
are  given  in  Figs.  68,  69,  and  JO  for  the  different  types  of  models  tested.  The  most 
accurate  set  of  assumptions  should  be  that  which  results  in  the  calculation  of 
Curve  c  in  Figs.  6j  through  70. 

A  comparison  of  measured  values  of  8  with  values  of  8  calculated  by  the  assump¬ 
tions  leading  to  Curve  c  is  given  in  Fig.  71*  The  scatter  of  data  in  this  figure  is 
very  large  relative  to  that  which  is  ordinarily  obtained  from  tests  of  a  flat  plate. 
For  instance,  the  root-mean-rquare  scatter  of  flat-plate  data  about  the  faired  line 
in  Fig.  15  is  only  3$  of  the  average  value  of  momentum  thickness  for  this  curve.  As 
noted  in  Fig.  71,  the  root -mean -square  scatter  of  data  is  approximately  11$  for  points 
half  way  along  the  curved  surface  and  approximately  17$  for  points  at  the  end  of  the 
curved  surface.  These  two  values  of  root -mean -square  scatter  would  be  reduced  to 
9$  and  13$  if  the  Mach  number  3*5  data  were  omitted.  The  average  value  of  the  ratio 
of  measured  8  to  calculated  8  is  0-99  at  a  point  half  way  along  the  curved  surface 
(upper  plot  of  Fig.  71)  and  is  1.04  at  the  downstream  end  of  the  curved  surface 
(lower  plot  of  Fig.  71 ) •  These  average  values  would  be  reduced  to  0.97  and  1.01  if 
the  Mach  number  3 ■ 5  data  were  omitted.  The  difference  between  the  measured  values 
and  a  value  of  unity  cannot  be  considered  significant  on  the  basis  of  the  large 
scatter  of  the  data  in  Fig.  71- 

Momentum  thicknesses  at  the  downstream  end  of  the  curved  surface  calculated  on 
the  basis  of  the  variation  of  G  with  length  determined  from  the  lag-length  procedure 
(Curve  e  in  Figs.  67  through  70)  are  different  from  those  calculated  on  the  basis  of 
the  measured  variation  of  G  with  length  (Curve  c).  The  average  value  of  the  ratio 
of  momentum  thickness  corresponding  to  Curve  e  to  that  corresponding  to  Curve  c  was 
O.98.  The  root -mean -square  difference  between  the  momentum  thicknesses  calculated 
by  the  two  different  procedures  was  4$  of  the  absolute  value  of  the  momentum  thick¬ 
ness  at  the  end  of  the  curved  portion  of  the  plate.  If  the  Mach  number  3*5  data 
were  omitted,  the  average  value  would  be  increased  to  0-99  and  the  root -mean -square 
scatter  could  be  reduced  to  3 $•  Therefore,  the  error  introduced  into  the  calculation 
of  growth  of  boundary  layer  momentum  thickness  by  the  fact  that  the  variation  of  G 
with  length  calculated  by  the  lag-length  procedure  is  different  from  the  actual 
variation  of  G  with  length  does  not  introduce  a  significant  error  into  the  calcuation 
of  boundary  layer  momentum  thickness. 
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All  the  data  presented  in  Figs.  67  through  71  were  obtained  with  zero  bleed.  On 
the  basis  of  theoretical  considerations,  boundary  layer  bleed  should  result  in  a  re¬ 
duction  in  boundary  layer  momentum  thickness.  If  the  flow  between  the  wall  and  a 
station  jq-,  away  from  the  wall  is  bled  off,  the  resulting  change  in  boundary  layer 
momentum  thickness  can  be  found  from  the  following: 

A6b"  ^v^yb^V(l"V)dy  (20) 

The  weight  flow  which  is  bled  off  is  given  by 

wb  =  yyb/3Vdy  (21) 


The  ratio  of  the  change  in  momentum  thickness  to  the  flow  which  is  bled  off  can  be 
obtained  through  dividing  Eq.  (20)  by  Eq.  (21); 


wb  yybyov  dy 


(22) 


Examination  of  this  equation  shows  that  the  greatest  change  in  boundary  layer  momen¬ 
tum  thickness  for  a  given  amount  of  bleed  occurs  if  the  flow  which  is  bled  off  has 
zero  velocity,  which  makes  the  term  within  the  parenthesis  in  the  preceding  equation 
equal  to  unity.  This  maximum  change  of  momentum  thickness  per  unit  bleed  flow  is, 
therefore,  given  as  follows: 


^tQ-THEORY 

wb 


(23) 


The  actual  change  of  boundary  layer  momentum  thickness  per  unit  bleed  flow  should  be 
less  than  that  given  by  the  preceding  equation  for  two  reasons:  (l)  the  flow  bled 
off  has  a  velocity  somewhat  greater  than  zero,  and  (2)  the  removal  of  the  low-energy 
flow  from  the  portion  of  the  boundary  layer  near  the  wall  should  result  in  an  increase 
in  wall  friction. 
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The  effect  of  boundary  layer  bleed  on  the  reduction  of.  boundary  layer  momentum 
thickness  was  determined  by  the  procedure  indicated  in  the  top  plot  of  Fig.  72.  The 
theoretical  variation  of  momentum  thickness  with  length  for  zero  bleed  given  in  this 
sketch  was  determined  by  the  procedure  for  Curve  c  in  Fig.  67 •  The  variation  of  G 
with  distance  for  this  theoretical  curve  was  obtained  by  fairing  between  a  value 
of  G  =  Gp  at  the  upstream  end  of  the  curved  surface  and  the  measured  value  of  G  at 
the  downstream  end  of  the  bleed  region.  It  was  not  possible  to  obtain  measured  values 
of  G  in  the  bleed  region  of  the  model  because  of  local  flow  angularity  caused  by  the 
flow  entering  the  perforations.  The  measured  difference  in  9  due  to  bleed, 

(  ^  0b  )  measured  was  as  the  difference  between  this  calculated  value  of  9 

and  the  measured  value  of  9  at  the  end  of  the  bleed  region. 

The  ratio  of  the  measured  decrease  in  9  per  unit  bleed  flow  divided  by  the 
theoretical  value  of  the  same  quantity  is  ,  iven  in  the  lower  plot  of  Fig.  J2 .  As 
noted  in  a  preceding  paragraph,  the  measured  decrease  in  9  due  to  bleed  should  be 
less  than  the  maximum  theoretical  value  of  this  quantity.  The  average  value  of  the 
ratio  of  the  measured  to  theoretical  decrease  in  9  due  to  bleed  is  on  the  order  of 
O.5..  The  relatively  large  scatter  of  the  small  amount  of  data  about  this  average 
value  results  from  inaccuracies  in  the  absolute  values  of  momentum  thickness  (see 
discussion  of  scatter  in  Fig.  7l). 


Boundary  Layer  Separation 

A  large  amount  of  data  on  the  pressure  ratio  for  incipient  separation  for  single¬ 
radius  two-dimensional  curved-surface  models  is  presented  in  Ref.  25.  These  data  and 
several  data  points  from  the  present  report  are  presented  in  Figs.  73  and  7^  as  a 
function  of  the  radius  of  curvature  of  the  model  divided  by  the  boundary  layer 
momentum  thickness  at  the  curved- surface  tangent  point.  Station  0.  The  arrows  on  the 
data  points  obtained  from  the  present  report  indicate  the  position  at  which  the 
points  should  fall  to  obtain  perfect  agreement,  with  the  data  from  Ref.  25.  The  agree¬ 
ment  between  the  data  from  the  present  report  and  from  Ref.  25  is  fair.  Also  shown 
in  these  figures  are  reference  curves  which  will  be  discussed  in  following  paragraphs. 
The  experimental  data  presented  in  Figs.  73  and  7^  indicate  that  the  pressure  ratio 
for  incipient  separation  increases  with  a  decrease  in  Reynolds  number  or  with  an  in¬ 
crease  in  W9q( decrease  in  dimensionless  pressure  gradient).  The  local  surface  angles 
and  Mach  numbers  corresponding  to  the  pressure  ratios  shown  on  the  left  side  of 
Figs.  73  and  7^-  are  given  cn  the  right  side  of  these  figures. 

In  order  to  explain  the  data  presented  in  Figs.  73  and  7^  in  terms  of  the  lag- 
length  concept  de/eloped  in  this  report,  a  series  of  calculations  were  performed 
with  systematic  variations  of  Mach  number,  initial  boundary  layer  momentum  thickness, 
and  Reynolds  number  based  on  initial  boundary  layer  momentum  thickness.  The  variation 
of  Mach  number  with  distance  along  the  curved  surface  was  assumed  to  be  that  deter- 
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mined  from  theory  (see  Fig.  66  for  the  effect  on  G  of  the  difference  between  ex¬ 
perimental  and  theoretical  Mach  number  distributions).  An  assumption  was  made  that 
separation  would  be  encountered  on  models  where  the  local  value  of  reduced  boundary 
layer  shape  parameter,  G,  exceeded  1.5.  This  value  of  G,  which  is  equal  to  a  value 
of  H  of  approximately  1.9  in  incompressible  flow,  was  selected  since  separation  i3 
usually  encountered  in  incompressible  flow  for  values  of  H  between  1.8  and  2.6.  The 
results  of  these  calculations  are  presented  as  dashed  lines  in  Figs.  73  and  7^*  The 
sharp  break  in  the  calculated  curves  near  values  of  W6q  of  150  for  initial  Reynolds 
numbers  of  2000  or  greater  occurs  because  the  lag  length  becomes  equal  to  the  length 
of  the  curved- surface  portion  of  the  model  at  these  break  points. 

The  effect  of  radius  of  curvature  and  Reynolds  number  on  the  pressure  ratio 
for  incipient  separation  calculated  from  the  lag-length  procedure  (the  dashed  lines) 
is  in  qualitative  agreement  with  the  data  obtained  from  Ref.  25.  The  reason  for 
the  differences  between  the  calculated  and  measured  curves  for  low  values  of  R /6q 
can  be  explained  on  the  basis  of  the  pressure  distributions  noted  in  Ref.  25  for  such 
models.  For  a  compression-comer  model (R/0O equal  to  zero),  the  increase  in  static 
pressure  due  to  the  comer  does  not  occur  in  zero  length  as  predicted  by  inviscid 
theory  but  spreads  out  over  a  considerable  distance  along  the  model.  The  pressure 
distribution  data  in  Ref.  25  show  that  the  actual  increase  in  static  pressure  on 
compression-corner  models  occurs  over  a  distance  which  is  considerably  greater  than 
the  lag  length  calculated  for  the  conditions  at  the  comer.  Therefore,  the  pressure 
ratio  for  incipient  separation  for  low  values  of  R/£?0  should  be  plotted  at  values  of 
R/0O  considerably  greater  than  zero.  Displacement  of  these  data  to  the  right  of 
the  break  points  on  Figs.  73  and  jk  would  result  in  considerable  improvement  in  the 
agreement  between  the  calculated  and  measured  pressure  ratios  for  incipient  separa¬ 
tion 


The  difference  between  the  measured  and  calculated  values  of  pressure  ratio 
for  incipient  separation  at  large  values  of  R /So  can  be  explained  on  the  basis  of 
the  theoretical  limiting  turning  angle  for  isentropic  external- compress ion  inlets. 
According  to  Ref.  40,  the  maximum  turning  angle  which  can  be  employed  in  external- 
compression  inlets  is  determined  by  limitations  on  the  shock  structure  where  the  waves 
from  the  compression  surface  coalesce.  This  limiting  turning  angle  is  noted  on 
Figs.  73  and  7^  IT  this  turning  angle  is  exceeded,  a  normal  shock  will  be  formed 
at  the  point  at  which  the  waves  from  the  compression  surface  coalesce.  A  normal 
shock  wave  was  noted  for  some  of  the  high-tuming-angle  models  on  which  data  were 
reported  in  Ref.  25;  therefore,  the  turning  angle  given  for  external-compression 
inlets  can  be  exceeded  to  some  degree  in  some  cases  without  causing  boundary  layer 
separation.  It  should  be  noted  that  pressure  ratios  considerably  greater  than  those 
determined  on  the  basis  of  the  limiting-turning-angle  criteria  could  be  obtained  if 
the  variation  of  static  pressure  with  distance  along  the  surface  were  created  in  a 
convergent  duct  instead  of  on  a  simple  curved  surface.  For  instance,  at  Mach  number 
3-5  the  pressure  ratio  for  incipient  separation  would  be  increased  from  12  for  the 
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simple  curved  surface  to  approximately  30  for  a  convergent  duct  at  large  values  of 
R /Qq  according  to  Fig.  7^* 

As  noted  in  a  preceding  section,  it  is  possible  with  a  given  model  to  obtain 
separation-free  flow  with  much  larger  values  of  boundary  layer  momentum  thickness  at 
Station  0  if  boundary  layer  bleed  is  employed  than  if  boundary  layer  bleed  is  not 
employed.  The  amount  of  boundary  layer  bleed  necessary  to  prevent  separation  on 
two  single-radius  two-dimensional  models  (R  =  1.5  and  3*0  in.)  is  given  in  the  upper 
plots  of  Figs.  75  and  j6  as  a  function  of  the  boundary  layer  momentum  thickness  at 
Station  0.  The  amount  of  bleed  necessary  to  prevent  separation  for  compression 
comers  (R  =  0),  calculated  on  the  basis  of  bleeding  all  of  the  boundary  layer  flow 
xrtiich  had  a  total  pressure  less  than  the  calculated  static  pressure  downstream  of  the 
compression  comer,  is  also  shown  in  the  upper  plots  of  Figs.  75  and  76. 

A  cross  plot  of  the  data  in  the  upper  plots  of  Figs.  75  and  76  is  given  on  the 
lower  plots  in  the  same  figures  as  a  function  of  dimensionless  gradient  (as  measured 
by  R/0O).  Although  it  was  not  possible  to  determine  any  significant  effect  of 
Reynolds  number  based  on  initial  momentum  thickness  for  the  Mach  2.5  data  shown  in 
Fig.  75,  there  appeared  to  be  a  trend  toward  increasing  the  bleed  flow  required  to 
prevent  separation  as  the  Reynolds  number  was  increased  at  Mach  3-0  (Fig.  76). 
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LIST  OF  SYMBOLS 


Cf 


G 

Gf 

Go 


Local  skin-friction  coefficient, 

(WALL  SHEAR  FORCE ) /( p^V2/  2 ) ( SURFACE  AREA) 

Reduced  boundary  layer  shape  parameter,  H/Hu-y 

Value  of  G  on  flat  plate,  1.12  (Rg/lOOO)'0'0^  (see  Fig.  l8) 

Value  of  G  at  upstream  end  of  adverse  pressure  gradient  in 
theoretical  calculations 


h  Tunnel  height  for  tests  of  two-dimensional  models  (see  Fig.  l) 

H  Boundary  layer  shape  parameter,  8 */ 8 

HN=7  i hn-5  Boundary  layer  shape  parameters  for  l/7  and  l/5-power  velocity 

profiles,  respectively  (see  Ref.  3l) 

J  Boundary  layer  deceleration  parameter  (see  Eq.  7) 

K  Match  ratio,  value  of  YmatCH / Yp  employed  in  calculation  of  Pg  (a 

value  of  K  =  1.3  was  employed  for  data  reduction  -  see  text) 

Lag  length  (see  Figs.  53  and  60) 

L|  Length  to  upstream  end  of  two-dimensional  curved- surface  model 

from  tunnel  throat  for  two-dimensional  models  (see  Fig.  l)  or 
from  model  leading  edge  for  axisyrametric  models  (see  Fig.  7) 

L  2  Length  from  upstream  end  of  curved- surface  model  to  Station  0 

(see  Figs.  3,  5,  and  7) 

L3  Circumferential  length  along  surface  of  radius  R  or  R^  (see 

Figs.  3,  5,  and  7) 

L4  Length  of  constant-angle  surface  between  two  curved  portions  of 

the  two-radii  two-dimensional  curved- surface  models  (see  Fig.  5) 

l_5  Circumferential  length  along  surface  of  radius  R2  (see  Fig.  5) 

L6  Length  along  constant-angle  surface  downstream  of  curved  portion 

of  model  (see  Figs.  3,  5,  and  j) 
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M 

M0 


Mg 

N 

n 

P 

pip 


MAX 


Sf 


>w 


MAX 


’W 


Mach  number  determined  from  Pgw/PijQ 

Nominal  tunnel  Mach  number  (also  Mach  number  at  upstream  end  of 
adverse  pressure  gradient  in  theoretical  calculations) 

Mach  number  at  outside  edge  of  boundary  layer 

1 

u  y  N 

Exponential  profile  parameter  in  equation  -jj-  =( -fc- ) 

Number  of  data  points  in  sample 
Pressure 

Injection  plenum  pressure 
Measured  pitot  pressure  (see  Fig.  9) 

Maximum  measured  pitot  pressure  (see  Fig.  9) 

Local  static  pressure  calculated  from  analytical  curve  fitted 
through  Pgw  and  Pgf  at  y^rpQjj  (see  text  and  Fig.  9) 

Static  pressure  calculated  from  measured  pitot  pressure,  Pp,  and 
tunnel  total  pressure,  P^q  (see  Fig.  9) 

Nominal  tunnel  static  pressure 

Maximum  static  pressure  measured  near  downstream  end  of  curved- 
surface  model 

Wall  static  pressure 

Local  total  pressure  calculated  from  Pp  and  Pg  (see  Fig.  9) 
Tunnel  total  pressure 


Dynamic  pressure. 


r 


Radius  of  axisymmetric  curved- surf ace  model  in  a  plane  normal 
to  axis  (see  Fig.  7) 
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re 


Maximum  radius  of  axisymmetric  curved- surface  model  in  a  plane 
normal  to  axis  (see  Fig.  7) 


Co  Radius  of  axisymmetric  curved- surface  model  at  Station  0  in  a 

plane  normal  to  axis  (see  Fig.  7) 

R  Radius  of  curvature  of  single-radius  models  (see  Figs.  3  and  7) 

R  Average  radius  of  curvature;  equals  R  for  single -radius  models 

and  (Rp  +  Rp)/2  for  two-radii  models 

R|  Initial  radius  of  curvature  of  two-radii  model  (see  Fig.  5) 


R2 

Rs 

R*o 


Tj 


u 

u' 

V 

vo 

vesp 

Wb 

WBL 


Second  radius  of  curvature  of  two-radii  model  (see  Fig.  5) 

Reynolds  number  based  on  boundary  layer  momentum  thickness. 

Values  of  Rg  at  Station  0  (also,  value  of  Rg  at  upstream  end 
of  adverse  pressure  gradient  in  theoretical  calculations) 

Total  temperature  -  assumed  constant  throughout  flow  field 
except  in  discussion  in  Appendix  I 

Local  velocity  determined  from  Pg/Pp 

Ideal  local  velocity  determined  from  Ps/PTq 

Velocity  determined  from  Pg^/Ppo 

Nominal  tunnel  velocity  (also,  velocity  at  upstream  end  of 
adverse  pressure  gradient  in  theoretical  calculations) 

Velocity  within  equivalent  step  profile  (see  Fig.  47) 

Flow  bled  from  the  boundary  layer 

/g 

pvdy 

o 

Flow  within  boundary  layer  at  Station  0 
Distance  from  tunnel  nozzle  throat 

Distance  from  tunnel  nozzle  throat  to  curved- surface  tangent 
point  (Station  0) 


H- 
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y 

Vb 

VmATCH 

yP 

z 

as 

aT 

°l 

a2 

/3 


y 

8 

So 
8* 
e  . 

90 

Adb 

(ADMEASURED 

(A0b)THEORY 


Distance  normal  to  surface  of  model 

/yb 

/jvdy 

Distance  normal  to  surface  where  y  =  K(yp)  (see  Fig.  9) 

Distance  normal  to  surface  at  which  Pp  is  a  maximum  (see  Fig.  9) 

Distance  across  model  (measured  from  centerline  of  tunnel) 

Local  surface  angle  (relative  to  surface  upstream  of  Station  0) 

Total  turning  angle  of  model  (see  Figs.  3,  5,  and  7) 

Initial  turning  angle  of  two-radii  model  (see  Fig.  5) 

Second  turning  angle  of  two-radii  model  (see  Fig.  5) 

Boundary  layer  flow  thickness,  defined  as  height  of  passage  which 
will  pass  flow  within  equivalent  step  profile  at  Mach  1.0  and 
free-stream  total  pressure  and  temperature  (see  Figs.  51  and  52 
and  Eq .  17 ) 

Ratio  of  specific  heats 

Boundary  layer  thickness,  determined  as  the  distance  normal  to 
the  surface  where  P>ji/PTq  = 

Boundary  layer  thickness  at  curved- surface  tangent  point  (Station  0) 
Boundary  layer  displacement  thickness  (see  Fig.  10) 

Boundary  layer  momentum  thickness  (see  Fig.  10) 

Boundary  layer  momentum  thickness  at  curved- surface  tangent  point 
(Station  0) 

Decrease  in  9  due  to  bleed  (see  Eq.  20) 

Measured  decrease  in  6  due  to  bleed  (see  Fig.  72) 

Theoretical  maximum  decrease  in  Q  due  to  bleed  (see  Fig.  72  and 

Eq-  23) 
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Local  density  determined  from  Pg/Pp 


Lag -length  residual,  [(f  )MEASURED‘  ( f  )  Eq  (I9)>F|G  J 

Density  determined  from  PSw/PTq 
Ideal  density  determined  from  Ps/Ptq 

Viscosity  evaluated  from  static  temperature  determined  from  V 
Bypass  angle  (see  Figs.  3  and  5) 

Conditions  at  beginning  of  integration  step  (see  Appendix  II) 
Conditions  at  end  of  integration  step  (see  Appendix  II ) 
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APPENDIX  I 


Effect  of  Wall  Temperature  on  Boundary  Layer  Characteristics 


The  literature  contains  a  large  amount  of  data  which  show  the  effect  of  wall 
temperature  on  boundary  layer  characteristics  in  flows  having  zero  pressure  gradient. 
For  instance,  in  Ref.  32  it  is  shown  that  a  decrease  in  wall  temperature  results  in 
an  increase  in  friction  coefficient  and  a  decrease  in  the  boundary  layer  shape 
parameter,  H,  relative  to  the  values  obtained  with  an  uncooled  wall.  Data  on  the 
effect  of  wall  temperature  on  friction  coefficient  and  reduced  boundary  layer  shape 
parameter  obtained  from  the  tables  in  Refs.  32  and  4l  are  given  in  Figs.  77  and  78* 

For  Mach  numbers  greater  than  zero,  the  values  of  G  obtained  from  Refs.  32  and  4l 
for  a  ratio  of  wall  temperature  to  insulated- wall  temperature  of  unity  are  less  than 
unity  (i.e.,  values  of  H  from  Refs.  32  and  4l  are  less  than  those  obtained  from 
Ref.  31  and  plotted  in  Fig.  13)*  This  difference  in  G  results  from  the  assumption 
of  a  recovery  factor  of  0-9  for  the  tables  in  Refs.  32  and  4l  as  compared  to  a  re¬ 
covery  factor  of  unity  for  the  tables  in  Ref.  31  • 

For  the  case  of  a  boundary  layer  flowing  over  a  cooled  wall,  the  energy  distri¬ 
bution  (total  temperature  distribution)  within  the  boundary  layer  is  determined  by 
the  rate  at  which  energy  is  transferred  to  the  wall  and  the  rate  at  which  energy  is 
transferred  to  the  boundary  layer  from  the  free  stream.  The  rate  of  transfer  of 
energy  from  the  free  stream  into  the  boundary  layer  is  much  greater  in  an  adverse 
pressure  gradient  than  in  zero  pressure  gradient  because  the  rate  of  transfer  of 
flow  from  the  free  stream  into  the  boundary  layer  is  much  greater  in  the  former  case 
than  in  the  latter.  Therefore,  the  effect  of  wall  temperature  on  friction  coefficient 
and  reduced  boundary  layer  shape  parameter  in  an  adverse  pressure  gradient  should  be 
much  less  than  that  shown  in  Figs.  77  and  jQ  for  zero  pressure  gradient. 

The  boundary  layer  momentum  equation  (Eq.  3)  applies  to  cooled  as  well  as  un- 
cogled  boundary  layers  if  actual  local  densities  are  employed  in  evaluating  the 
boundary  layer  momentum  and  displacement  thicknesses.  If  the  use  of  a  cooled  wall 
has  anywhere  near  as  large  an  effect  on  G  in  an  adverse  pressure  gradient  as  in  zero 
pressure  gradient  (see  Fig.  j8) ,  then  the  boundary  layer  growth  in  an  adverse  pressure 
gradient  with  cooling  should,  according  to  Figs.  42,  43,  and  44  be  much  less  than  the 
boundary  layer  growth  in  an  adverse  pressure  gradient  without  cooling.  Therefore, 
determination  of  the  effect  of  wall  cooling  on  boundary  layer  characteristics  in  an 
adverse  pressure  gradient  is  extremely  important. 

As  noted  in  the  preceding  paragraph,  the  effect  of  wall  temperature  on  boundary 
layer  characteristics  is  different  in  an  adverse  pressure  gradient  than  in  zero 
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pressure  gradient  because  of  the  large  influx  of  flow  into  the  boundary  layer  in  an 
adverse  pressure  gradient.  One  method  of  analyzing  the  problem  of  the  determination 
of  the  effect  of  wall  cooling  on  boundary  layer  characteristics  in  an  adverse 
pressure  gradient  might  be  to  adopt  a  step  profile  in  which  the  total  temperature  of 
the  step  profile  was  a  variable  rather  than  a  value  assumed  equal  to  free- stream 
total  temperature.  The  total  temperature  within  the  equivalent  step  profile  could 
then  be  determined  by  an  energy  balance  at  each  station  along  the  surface  in  the 
adverse  pressure  gradient. 
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APPENDIX  II 

Development  of  Incremental  Form  of  Boundary  Layer  Growth  Equation 

The  differential  form  of  the  houndary  layer  growth  equation,  Ekj.  (12),  is  as 
follows : 

d0  =  tj-  dx-SJdM  (24) 

Rearranging  Ekj.  (24)  yields 


de=(|j  m~e) JdM  (25) 


This  equation  is  to  he  integrated  over  a  finite  change  in  Mach  number.  During  this 
integration, cf  ,  dx/dM  and  J  are  assumed  to  be  constant.  Then 
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02-£i  dx 
2  2J  5m 


— - |--j~  =-J(M2-M,)=-JAM 

9{  ~2J  dM 


e  q  dx 
2~  2 J  dM  ..JAM 

01_£f  f!*'e 


2 J  dM 


Expanding  the  right  side  of  Eq.  (29)  in  a  power  series  yields: 


a  *1 

2_2J  dM 


0,-21 

1  2  J 


dx 

dM 


*i  + 


(-JAM)  [l 


,  (-JAM)  (-JAM)  , 
2  6  + 


„  Cf  dx  n  Cf  dx  ( cf  dx  ,  ,  „ 

e2  -  2J  dM  =  ®I  ■  iz  Bm+  (  2T  BM  (' JAM)_  B'  JAM 


r  JAM  ,  (JAM)"  ,  1 

L l~~s~  +  —  +  •  ‘  ‘J 


(28) 


(29) 


(30) 


(31) 
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Since  dx/dM  has  been  assumed  constant  during  a  step. 


d  x  .  Ax 

"3m  '  SW 


Equation  (32)  then  becomes 


A  Q  -  ^  A  x  -  £?|  JAM  )  |)  ■ 


JAM  ,  (JAM)2. 
2  *  6  + 


Equation  (3^)  is  similar  to  Eq.  (12)  except  for  the  second  term  on  the  right  side  of 
the  equation,  which  is  a  terra  to  be'  used  when  A 8  is  an  appreciable  fraction  of  8\ . 
For  most  calculations,  Eq.  (13),  which  is  a  reduced  form  of  Eq.  (3*0,  contains 
enough  terms . 


A 0=  (  -^AX  —  ^  JAM) 


R-I285-  1 1  TABLE  I 
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TABLE  H 


COMPARISON  OF  MEASURED  VALUES  OF  LAG  LENGTH 
WITH  VALUES  COMPUTED  FROM  FITTED  CURVE 


Pg  -  (4/9)  MEASURED  —  (i/©  )[EQ.(  19  ),  FIG.60J 


TYPE  OF  DATA 

NUMBER  OF 
POINTS,  n 

AVERAGE  LAG -LENGTH 
RESIDUAL,  Ifig/i 1 

ALL 

1  49 

0 

SINGLE-RADIUS 
TWO-DIMENSIONAL  MODELS 

82 

-2.9 

TWO-RADII 

TWO-DIMENSIONAL  MODELS 

44 

+  4.0 

AXISYMMETRIC  MODELS 

23 

+  2.7 

R/e0  FROM  0  TO  250 

40 

+  3.0 

n/e0  FROM  2  50  —  300 

39 

-3.8 

r/B0  FROM  300  —  400 

36 

-0.8 

R /Bq  FROM  400  —  700 

34 

+  1 .8 

M0  =  2.0 

1  1 

-3.3 

M  o  =2,5 

52 

-  1 .4 

MO  =  3.0 

:  " 

4  3.9 

Mq  *  3.5 

12 

-1  5.3 

46 


VALVt  FOR  CONTROLLING  PRESSURE 
IN  INJECTION  PLENUM 


FLOW  INJECTION  - 
REGULATING  VALVE 


PHOTOGRAPH  OF  SINGLE-RADIUS  TWO-DIMENSIONAL  CURVED -SURFACE  MODELS 
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FIG. 6 
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AXISYMMETRIC  CURVED-SURFACE  MODELS 


PHOTOGRAPH  OF  AXISYMMETRIC  CURVED  -  SURFACE  MODELS 
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FIG  8 
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STATIC  PRESSURE  Ps  PRESSURE 

TUNNEL  TOTAL  PRESSURE’  PT„  TUNNEL  TOTAL  PRESSURE 
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FIG. 10 


EXAMPLE  OF  MASS  FLOW  AND  MOMENTUM  FUNCTIONS 
FOR  BOUNDARY  LAYER  HAVING  LARGE  STATIC 
PRESSURE  GRADIENT  NORMAL  TO  WALL 
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TOTAL  TEMPERATURE  ASSUMED  CONSTANT  THROUGHOUT  BOUNDARY  LAYER 


■P  U,  IDEAL  MASS  FLOW  PER  UNIT  AREA 
CALCULATED  FROM  PC/PT 


P  V,  IDEAL  MASS  FLOW  PER  UNIT  AREA 
-  W 

CALCULATED  FROM  Pc  /P_ 


MATCH 


P U,  ACTUAL  MASS  FLOW  PER  UNIT  AREA 
CALCULATED  FROM  Pc  /  P_ 


< 

LU 

cr 

< 


cc 

LU 

a. 


o 


p'u'2,  IDEAL  MOMENTUM  PER  UNIT  AREA 


CALCULATED  FROM  Pc 


^WV  ,  IDEAL  MOMENTUM  PER  UNIT  AREA 
CALCULATED  FROM  Ps  /  PT 


full',  ACTUAL  MASS  FLOW  PER  UNIT  AREA 
CALCULATED  FROM  PSL/PPL  AND 
MULTIPLIED  BY  IDEAL  VELOCITY 
CALCULATED  FROM  Pc  /P-c 


y match 


Fu  ,  ACTUAL  MOMENTUM  PER  UNIT  AREA 
CALCULATED  FROM  Ps/Pp 


DISTANCE  NORMAL  TO  SURFACE  OF  MODEL,  y  -  IN 
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WALL  STATIC  PRESSURE 


TYPICAL  WALL  STATIC  PRESSURE  DISTRIBUTIONS 
ON  REFERENCE  FLAT  PLATE 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS  12  AND  14 


STATIC  PRESSURE  DISTRIBUTION!  ACROSS  PLATE,  X=I5.6  IN 


RIGHT  SIDE  WALL 
LOOKING  DOWNSTREAM 


LEFT  SIDE  WALL 

LOOKING  DOWNSTREAM  DISTANCE  ACROSS  MODEL,  Z  -  IN 


CENTERLINE  STATIC  PRESSURE  DISTRIBUTION 
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DISTANCE  FROM  TUNNEL  NOZZLE  THROAT,  X-IN 


BOUNDARY  LAYER  SHAPE  PARAMETER  FOR  l/ 7  —  POWER  PROFILE, 
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BOUNDARY  LAYER  SHAPE  PARAMETER  FOR  l/7-P0WER  PROFILE 

TOTAL  TEMPERATURE  ASSUMED  CONSTANT  THROUGHOUT  BOUNDARY  LAYER 
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FIG  14 


TYPICAL  VALUES  OF  REDUCED  BOUNDARY  LAYER 
SHAPE  PARAMETER  FOR  REFERENCE  FLAT  PLATE 

M0=  3.0  PIP/PT0;0I° 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  II  AND  12 
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BOUNDARY  LAYER  MOMENTUM  BOUNDARY  lAYER  DISPLACEMENT  BOUNDARY  LAYER  THICKNESS 

THICKNESS,  0  -  IN  THICKNESS,  8*  -IN 


EFFECT  OF  AIR  INJECTION  ON  BOUNDARY  LAYER 
THICKNESS  FOR  REFERENCE  FLAT  PLATE 


INJECTION  PLENUM  PRESSURE 


TUNNEL  TOTAL  PRESSURE 
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FIG  17 


EFFECT  OF  AIR  INJECTION  ON  BOUNDARY  LAYER 
VELOCITY  PROFILES  FOR  REFERENCE  FLAT  PLATE 
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FIG.  19 


TYPICAL  SCHLIEREN  PHOTOGRAPHS  OF  SINGLE  -  RADIUS 
TWO-DIMENSIONAL  CURVED  -  SURFACE  MODEL 

Mq  >  3  R  :  6  IN.  a  T  :  30  DEG  Wb  =  0 


V 

>  STREAKS  IN  TUNNEL 

WINDOW  GLASS 

^VARIABLE  DIFFUSFR 

NOZZLE  BuOCK 

f 

■  . 

, 

- 

✓  ' 

TANGENT  POINT  * 

STAT ION  0 

TEST  MODEL 

L  j  3  1  A  IN 

L  6  -  0.20  IN 

®0 

=  0.0096  IN.  RgQ=l530 

80=O.I7  in. 

t 

v» 

dQ  >0.0205  IN.  3250  SQ  >0.33  IN. 
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TYPICAL  WALL  STATIC  PRESSURE  DISTRIBUTIONS 
ON  SINGLE -RADIUS  TWO-DIMENSIONAL 
CURVED -SURFACE  MODEL 


Mo  =  3’°  R  =  6  IN.  <2  T  =  30  DEG  =0 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  21, 22 , 23,24,46,53, 55,62,  AND  68 
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FIG  21 
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LOCAL  SURFACE  ANGLE,  Qj— DEG 


TYPICAL  VARIATION  OF  MACH  NUMBER  WITH  DISTANCE  ALONG 
SINGLE-RADI'JS  TWO-DIMENSIONAL  CURVED- SURFACE  MODEL 


R  =  6  IN. 


Ot  -  30  DEG 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS  2  0, 22 , 23 , 24,4  6, 53,  55 , 62 ,  AND  68 


M  IS  MACH  NUMBER  DETERMINED  FROM  Pc  /PT 

SW/  T0 

Mg  IS  MACH  NUMBER  AT  OUTSIDE  EDGE  OF  BOUNDARY  LAYER 


SYMBOL 

z 

o 

Q> 

80-'N 

Rfl 

°0 

>-o  ' 

o 

A 

- - 

0  0105 

0  0160 

'0  02  i  5 

0.17 

0  30 

0  33 

1680 

2540 

3420 

30 


^  28 

a 

°  2  6 

S' 


2  2 


2  O 


?  -I  0  I  2  3 

DISTANCE  FROM  CURVED- SURFACE  TANGENT  POINT,  (X-XQ)-IN. 


67 


DISTANCE  NORMAL  TO  SURFACE  OF  MODEL 
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FIG. 22 


TYPICAL  BOUNDARY  LAYER  VELOCITY  PROFILES  FOR 
SINGLE- RADIUS  TWO-DIMENSIONAL  CURVED-SURFACE  MODEL 

M0=3  0  R  =  6  IN  aT  =  30  DEG  wb=  0 

0-  =  0.0160  IN.  R a  ■  2540 

O  °0 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  20,21,23,24,46,53,55,62,  AND  68 
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FIG  23 


EFFECT  OF  APPROACH  BOUNDARY  LAYER  THICKNESS  ON  VELOCITY 
PROFILES  AT  DOWNSTREAM  END  OF  ADVERSE  PRESSURE  GRADIENT 
FOR  SINGLE -RADIUS  TWO-DIMENSIONAL  CURVED -SURFACE  MODEL 

M0=  3.0  R  =  6  IN.  as=aT  =30  OEG 

Wb  =  0  X  -  XQ  =  3.14  IN. 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  20,  21 ,22,24,46,53,55,  62,  AND  68 
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CONDITIONS  AT  ds  =  30  DEG 
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^-BOUNDARY  LAYER  SEPARATED  NEAR  ds  =  15  DEG 
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DISTANCE  FROM  CURVED- SURFACE  TANGENT  POINT,  (x  -  X0) — IN 


LOCAL  SURFACE  ANGLE,  us-DEG 


TWO-DIMENSIONAL  CURVED-SURFACE  MODEL 


Mq  =  30  R  =  6  IN  aT  =  30  DEG  Wb  =0 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  20,21, 22,23,46,53,55, 62,  AND  68 


TYPICAL  BOUNDARY  LAYER  CHARACTERISTICS  FOR  SINGLE- RADIUS 
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(  INCIPIENT  SEPARATION  ) 
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FIG.  27 


TYPICAL  SCHLIEREN  PHOTOGRAPH  OF  AXISYMMETRIC 
CURVED  -  SURFACE  MODEL 


Mo  =  2.5  R  =  3  IN.  ar=30  DEG  Wb  =  O 

6q  -  0.0072  IN.  R0O=I5IO  S0  =  0.11  IN. 
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TYPICAL  WALL  STATIC  PRESSURE  AND  MACH  NUMBER 
DISTRIBUTIONS  ON  AXISYMMETRIC  CURVED  -  SURFACE  MODEL 

Mo  =  3.0  .  R--I.51N.  aT  ~  30  DEG  Wb=  0 

°  80  =0.0056  IN.  R0Or900 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  29,59,64,  AND  70 
DIFFERENT  SYMBOLS  INDICATE  DATA  OBTAINED  ON  DIFFERENT  DAYS 


MEASURED 


ME0R6 


DISTANCE  FROM  CUR  VE  D  -  SURFACE  TANGENT  POINT,  (  X-X0)  -  IN 


LOCAL  SURFACE  ANGLE, a 
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REDUCED  GOUNOARY  LAYER  SHAPE  BOUNDARY  LAYER  MOMENTUM  BOUNDARY  LAYER  DISPLACEMENT 
PARAMETER,  G  -  H/Hn,7  THICKNESS,  0-IN  THICKNESS,  5*- IN 
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FIG.  29 


TYPICAL  BOUNDARY  LAYER  CHARACTERISTICS 
FOR  AXISYMMETRIC  CURVED -SURFACE  MODEL 

-  3.0  R  =  I  5  !N.  a  T=  30  DEG  W^O 

0O;O-OO56IN.  R0or  900 

ADDITIONAL  OATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  28  ,  59 , 64,  AND  70 
DIFFERENT  SYMBOLS  INDICATE  DATA  OBTAINED  ON  DIFFERENT  DAYS 
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FIG.  30 
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C)  INCIPIENT  SEPARATION  ( Wb/*0Lo  s°  lS®  1 


R-I285-II 


FIG.  31 


SCHLIEREN  PHOTOGRAPH  OF  TWO-DIMENSIONAL  CURVED-SURFACE 
MODEL  WITH  MORE  BOUNDARY  LAYER  BLEED  THAN 
THAT  REQUIRED  TO  PREVENT  SEPARATION 


Mr.5  2.5 


R  =  3  IN. 


a  T  =  30  DEG 


8Q  =0.0146  IN. 


Rg  =  3050 
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MACH  NUMBER,  M  CR  Ms  WALL  STATIC  PRESSURE 
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FIG.  32 


TYPICAL  WALL  STATIC  PRESSURE  AND  MACH  NUMBER 
DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  CURVED-SURFACE 
MODEL  WITH  BOUNDARY  LAYER  BLEED 


M 
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THEORY 
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MEASlljREO  M 


ORtflCALllWVj^ 


IcjuNUARYT  uEHSl 


3  0  R  =  3.0  IN.  aT  =  30  DEG  0O  =  O.O228IN.  R§0  =5620 

Wb  /  WBL0  =  O134  (Wb/WBL0  :°  078  AT  INCIPIENT  SEPARATION) 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  33  AND  34 
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DISTANCE  FROM  CURVED- SURFACE  TANGENT  POINT,  (x -XQ)  -  IN. 
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LOCAL  SURFACE  ANGLE,  as~  DEG 
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DISTANCE  NORMAL  TO  SURFACE  OF  MODEL, 
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FIG.  33 


TYPICAL  BOUNDARY  LAYER  VELOCITY  PROFILES  FOR 
TWO-DIMENSIONAL  CURVED- SURFACE  MODEL 
WITH  BOUNDARY  LAYER  BLEED 

Mq  =  3.0  R=  3.0  IN.  Of  =  30  DEG  Q0=0.0228IN.  -  3620 

Wb/ WBl0  =  °-134  (  W  b/WB1_0  =  0.078  AT  INCIPIENT  SEPARATION) 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS  32  AND  34 
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REDUCED  BOUNDARY  LAYER  SHAPE  BOUNDARY  LAYER  DISPLACEMENT 

PARAMETER,  G =  h/hn=7  BOUNDARY  LAYER  MOMENTUM  THICKNESS,  8*- IN. 

THICKNESS,  0  -  IN 


R - 1285  - 1 1 


FIG  34 


TYPICAL  BOUNDARY  LAYER  CHARACTERISTICS  FOR 
TWO-DIMENSIONAL  CURVED  -  SURFACE  MODEL 
WITH  BOUNDARY  LAYER  BLEED 

M0  -  3,0  f):  3.0  IN.  ClT=30DEG  #O  =  OO220IN-  R0O  =  362° 

wb/  WBL0  =0.134  (wb/WB|_0  =0.078  AT  INCIPIENT  SEPARATION  ) 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  32  AND  33 


DISTANCE  FROM  CU R V E D  -  SUR FAC E  TANGENT  POINT,  (  X-X0  )  -  IN 
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DISTANCE  NORMAL  TO  SURFACE  OF  MODEL, 
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RG.35 
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FIG.  36 
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200  500  1000  2000  5000  10,000  20,000  50,000 

REYNOLDS  NUMBER  BASED  ON  BOUNDARY  LAYER  MOMENTUM  THICKNESS , Rg 


LOCAL  SKIN -FRICTION  COEFFICIENT 
LOCAL  SKIN  -  FRICTION!  COEFFICIENT  FOR 
FLAT  PLATE  AT  SAME  REYNOLDS  NUMBER 


R  -  1285  -  I  I 


FIG  37 


EFFECT  OF  REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER 
ON  LOCAL  SKIN- FRICTION  COEFFICIENT 


CURVES  OBTAINED  FROM  REFS.  26  AND  3? 


(SEE  TEXT) 
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REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER  G 

REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER  ’  GF 
FOR  FLAT  PLATE  AT  SAME  REYNOLDS  NUMBER 


83 


R- 1285  -  1 1 


FIG.  39 


BOUNDARY  LAYER  DECELERATION  PARAMETER 
FOR  MACH  NUMBERS  0.5  TO  18 


j  = 


dq/q  .  rH  d  v/vi 

d  m  N  - 7  d  M  J 


J  USED  IN  FOLLOWING  EQUATION  (EQ.  12  ) 
A  9  -■  %}-  AX  -  0JAM 
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BOUNDARY  LAYER  DECELERATION  PARAMETER 
FOR  MACH  NUMBERS  1.5  TO  3.5 


FIG  40 
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FIG  4  <2 


VARIATION  OF  BOUNDARY  LAYER  MOMENTUM  THICKNESS  WITH 
MACH  NUMBER  FOR  INITIAL  MACH  NUMBER  3.0 


ZERO  //ALL  FRICTION  ASSUMED 


MACH  NUMBER,  M 


8  8 
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FIG  43 


VARIATION  OF  BOUNDARY  LAYER  MOMENTUM  THICKNESS  WITH 
MACH  NUMBER  FOR  INITIAL  MACH  NUMBER  5.0 


ZEf<0  WALL  FRICTION 


ASSUMED 


MACH  NUMBER,  M 
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BOUNDARY  LAYER  MOMENTUM  THICKNESS _  8 

BOUNDARY  LAYER  MOMENTUM  THICKNESS  FOR  M  =  10  ’  0M  =  io 
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FIG  4  4 


VARIATION  OF  BOUNDARY  LAYER  MOMENTUM  THICKNESS  WITH 
MACH  NUMBER  FOR  INITIAL  MACH  NUMBER  10.0 


ZERO  WALL  FRICTION  ASSUMED 


NUMBER 
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FIG  45 


VARIATION  OF  TOTAL  FLOW  WITHIN  BOUNDARY  LAYER 
WITH  MACH  NUMBER  FOR  INITIAL  MACH 
NUMBERS  3.0 ,5.0,  AND  10.0 


G  =  10 

ZERO  WALL  FRICTION  ASSUMED 


CD 
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a: 
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MACH  NUMBER,  M 
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FLOW  WITHIN  BOUNDARY  LAYER  AT  STATION 


R- 1285- 1 1 


FI  6. 47 


COMPARISON  OF  EQUIVALENT  STEP 
ACTUAL  BOUNDARY  LAYER 


PROFILE 

PROFILE 


AND 


ACTUAL  BOUNDARY 
LAYER  PROFILE 


EQUIVALENT  STEP 
PROFILE 


NOTES : 

EQUIVALENT  STEP  PROFILE  AND  ACTUAL  BOUNDARY  LAYER  PROFILE 
HAVE  SAME  VALUE  OF  8*  AND  9  AND,  HENCE,  H  AND  G. 

BOUNDARY  LAYER  THICKNESS,  8,  AND  FLOW  WITHIN  BOUNDARY 

LAYER,  /SpVdY,  ARE  LESS  FOR  EQUIVALENT  STEP 
o  r 

PROFILE  THAN  FOR  ACTUAL  BOUNDARY  LAYER  PROFILE- 

LOCAL  STATIC  PRESSURE  IS  INDEPENDENT  OF  DISTANCE  NORMAL 
TO  WALL  FOR  EQUIVALENT  STEP  PROFILE. 
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REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER 


R-  1285-  I  I 


FIG  48 


EFFECT  OF  VELOCITY  WITHIN  EQUIVALENT  STEP  PROFILE 
ON  REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER 


R-I285-  II 


FIG  50 


EFFECT  OF  RAPID  DECREASE  IN  MACH  NUMBER  ON 
REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER 
FOR  COMPRESSIBLE  FLOW 

Y  -  I.  4 


ASSUMPTIONS 

1.  CONSTANT  TOTAL  PRESSURE  WITHIN 
STEP  PROFILE 

2.  ZERO  MIXING  FROM  OUTSIDE 
BOUNDARY  LAYER 

3.  ZERO  WALL  FRICTION 
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FIG.  51 


r_  1285 -  1 1 

VARIATION  OF  DIMENSIONLESS  FLOW  THICKNESS  WITH 


REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER 
FOR  MACH  NUMBERS  0.6  TO  4.5 


boundary  layer  flow  thickness 
bouncary  layer  momentum  thickness  ’ 


TYPICAL  EXAMPLE  OF  VARIATIONS  OF  ACTUAL  AND  IDEAL  BOUNDARY 
LAYER  FLOW  THICKNESS  AND  LAG  LENGTH  WITH  DISTANCE  ALONG 
SINGLE  -  RADIUS  TWO-DIMENSIONAL  CURVED -SURFACE  MODEL 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND 
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FIG  54 


SIMPLIFIED  METHOD  FOR  DETERMINATION  OF  REDUCED 
BOUNDARY  LAYER  SHAPE  PARAMETER  FOR  SMALL 
ADVERSE  MACH  NUMBER  GRADIENTS 
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DETERMINATION  OF  EFFECTIVE  VALUE  OF  M0-  M 


DISTANCE 
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c/1 


DETERMINATION  OF  G  FROM  EFFECTIVE  VALUE  OF  M0-M 
(SEE  FIG.  50) 
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R-I285-II 


FIG. 55 
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FIG  56 


TYPICAL  VALUES  OF  LAG  LENGTH  FOR  SINGLE-RADIUS 
TWO-DIMENSIONAL  CURVED- SURFACE  MODEL 

M0  =  3,0  R  =  6  IN.  a j  =  30  DEG  wb  =  0 
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REYNOLDS  NUM8ER  BASED  ON  BOUNDARY  LAYER  MOMENTUM  THICKNESS,  Rg 
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FIG  57 


R-  1285  -  1 1 

TYPICAL  VALUES  OF  BOUNDARY  LAYER  FLOW  THICKNESS 
FOR  TWO-RADII  TWO-DIMENSIONAL  CURVED  -  SURFACE  MODEL 

Mq  =  3.0  R  |  =  S  IN.  R2  =  3IN.  a  |  =  20  DEG  a  T  =35  OEG  Wb:° 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  F I6S  26,  63,  AN  D  fi9 
ARROWS(t  )  INDICATE  STATIONS  AT  WHICH  VALUES  OF  -^/f?GIVEN  IN  FIG  b8  WERE  DETERMINED 
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FIG  58 


TYPICAL  VALUES  OF  LAG  LENGTH  FOR  TWO- RADII 
TWO-DIMENSIONAL  CURVED -SURFACE  MODEL 


M0=  3.0 


R I  =  6  IN. 


R  2  =  3  IN. 


a | =20  DEG 


aT=35  DEG 


Wb=  0 
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BOUNDARY  LAYER  FLOW  THICKNESS 


R-I2B5-II 


FIG.  59 


TYPICAL  VALUES  OF  BOUNDARY  LAYER  FLOW  THICKNESS  AND 
LAG  LENGTH  FOR  AXISYMMETRIC  CURVED-SURFACE  MODEL 


M0=3.0 


R  1.5  IN 
Q  Q  -  0.0056  IN. 


a  T  =  3  0  DEG 


Ro  =-900 
fo 


wb=  0 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  28 , 29.  64 ,  AND  70 
ARROWS  U)  INDICATE  STATIONS  AT  WHICH  VALUES  OF  £/B  WERE  DETERMINED 
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LAG  —  LENG  1  H  RESIDUAL  ,  p  z  /  0  wE  A  SORED  ”  ^ 
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FIG  61 


COMPARISON  OF  LAG  LENGTHS  FROM  DATA 
AND  FROM  FITTED  CURVES 

AVERAGE  VALUE  =  0 

ROOT  -  MEAN  -  SQUARE  SCATTER  =  15.8 
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REYNOLDS  NUMBER  BASED  ON  BOUNDARY  LAYER  MOMENTUM  THICKNESS , Rg 
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FIG  62 


TYPICAL  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES 
OF  REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER  FOR 
SINGLE-RADIUS  TWO-DIMENSIONAL  CURVED-SURFACE  MODEL 

M0=  3.°  R-6  IN.  CrTs  30  DEG 


EASU 


jCALpOLAte] 


4C.WATI3D 


ALjCUL 


DISTANCE  FROM  CURVED- SURFACE  TANGENT  POINT,  (X-XQ) -IN 


— -  0  5  10  15  20 

LOCAL  SURFACE  ANGLE,  <2=—  DEG 


6  or- 

0.0 

160  IN 

| 

*ff0 

T 

*25*0 

"TT 

1 

f 

; 

CAliCUl 

ATED  ^ 

METHOI 

OF 

-4 

r 

. 

■ 

*  i 

L 

-  — - 

en=!o.o«o9 

4" 

RoJs|6hO 

_ 

'CALCUL 

ATEJD  6 

Y 

METHOC 

_ L. 

OF 

A€F 

26- 

k 

r  — 

A 

k 

^  H-f 

V 

7 

1 

<21 

(X- 

V 

>  1.35 

r  t  r 

i 

q 

. 

' _ 

— . 

RON 

v; 

“"oi 

^  FtO  610 

108 


TYPICAL  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF 
REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER  FOR  TWO- RADII 
TWO-DIMENSIONAL  CURVED -SURFACE  MODEL 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS  26,  57, AND  69 


DISTANCE  FROM  CURVED-SURFACE  TANGENT  POINT, 


LOCAL  SURFACE  ANGLE,  CL 


R-I285-II 


FIG  64 


TYPICAL  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES 
OF  REDUCED  BOUNDARY  LAYER  SHAPE  PARAMETER  FOR 
AXI SYMMETRIC  CURVED -SURFACE  MODEL 


M0  =  3.0 


So 


R  =  1.5  IN. 

0  0056  IN. 


a  T  =  30  DEG 

R0O  s  900 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  28,29,59,  AND  70 
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FIG. 65 


SUMMARY  OF  DATA  ON  DIFFERENCE  BETWEEN 
CALCULATED  AND  MEASURED  MAXIMUM  VALUES 
OF  REDUCED  BOUNDARY  LAYER 
SHAPE  PARAMETER 
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FIG.  66 


EFFECT  OF  SMALL  CHANGES  IN  MACH  NUMBER  DISTRIBUTION 
ON  CALCULATED  VALUES  OF  REDUCED  BOUNDARY  LAYER 
SHAPE  PARAMETER  FOR  SINGLE  -  RADIUS 
TWO-DIMENSIONAL  CURVED-SURFACE  MODEL 
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FIG.  67 


ASSUMPTIONS  EMPLOYED  IN  CALCULATING  BOUNDARY  LAYER 
MOMENTUM  THICKNESS  GROWTH  ON  CURVED  -  SURFACE  MODELS 


CURVE 

IDENTIFICATION 

G  USED  IN 
CALCULATION 

Cf  USED  IN 
CALCULATION 

a 

Gp  ,  FIG.  18 

0 

b _ _ 

MEASURED 

0 

c 

MEASURED 

FROM  FIGS.  36 
AND  37 

d 

MEASURED 

FROM  FIG.  36 

e 

CALCULATED  FROM 
X/B  IN  FIG.  60 

FROM  FIGS.  36 
AND  37 

o 

DISTANCE  FROM  CURVED  -  SURFACE  TANGENT  POINT, (X-X0) 


8b~9Q  CAUSED  BY  FACT  THAT  ACTUAL  G  IN  ADVERSE  PRESSURE  GRADIENT 
IS  GREATER  THAN  G  THAT  WOULD  EXIST  ON  FLAT  PLATE  AT  SAME  R0 

8c~8b  CAUSED  BY  WALL  FRICTION 

8d-8c  CAUSED  BY  FACT  THAT  FRICTION  IN  ADVERSE  PRESSURE  GRADIENT  IS 
LESS  THAN  FRICTION  ON  FLAT  PLATE  BECAUSE  OF  REDUCED  VELOCITY 
GRADIENT  AT  WALL 

9e~8c  CAUSED  BY  FACT  THAT  VALUES  OF  G  CALCULATED  BY  LAG-LENGTH 
PROCEDURE  ARE  NOT  EXACTLY  EQUAL  TO  MEASURED  VALUES  OF  G 


I  I  3 


R  -1285  -I  I 


FI  G.  68 


OISTANCE  FROM  CURVED  -  SURFACE  TANGENT  POINT,  (X-X0)-IN 
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TYPICAL  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF 
BOUNDARY  LAYER  MOMENTUM  THICKNESS  FOR  SINGLE-RADIUS 
TWO-DIMENSIONAL  CURVED- SURFACE  MODEL 


M0=  3.0  R  =  6  IN.  aT  =  30  DEG  Wb  =  0 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  20, 21,22, 23 , 24 , 4  6, 53, 55,  AND  62 

-  MEASURED 

}  CALCULATED  (SEE  FIG  67) 
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FIG  69 


TYPICAL  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF 
BOUNDARY  LAYER  MOMENTUM  THICKNESS  FOR  TWO-RADII 
TWO-DIMENSIONAL  CURVED-SURFACE  MODEL 

MQ=3.0  R  |  =  6  !N.  R2  =  3  1N.  a,  =20  DEG  a  T  =35  DEG  W|j=0 

ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  26,57,  AND  63 
-  MEASURED 

|  CALCULATED  (SEE  FIG.  67) 
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FIG.  70 


TYPICAL  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF 
BOUNDARY  LAYER  MOMENTUM  THICKNESS  FOR  AXISYMMETRIC 

CURVED -SURFACE  MODEL 

M0=  3.0  R  =  1.5  IN.  aT  30  DEG  =  ° 


MEASURED 

CALCULATED  (SEE  FIG  67) 


ADDITIONAL  DATA  FOR  THIS  MODEL  FOUND  IN  FIGS.  28,  29,  59,  AND  64 
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FIG  71 


COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES 
OF  BOUNDARY  LAYER  MOMENTUM  THICKNESS 
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FIG  72 


EFFECT  OF  BOUNDARY  LAYER  BLEED  ON  GROWTH  OF 
BOUNDARY  LAYER  MOMENTUM  THICKNESS 


DISTANCE  FROM  CURVE D  -  SURFACE  TANGENT  POINT,  X -X0 
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PRESSURE  RATIO  FOR  INCIPIENT  SEPARATION  ON  SINGLE -RADIUS 
TWO-DIMENSIONAL  CURVED -SURFACE  MODEL 
AT  MACH  NUMBERS  2.0  AND  2.5 


wb=o 

CALCULATED  FROM  4/8  DATA  IN  FIG.  60  ASSUMING  SEPARATION  OCCURS  AT  G=  1.5 
MAXIMUM  VALUE  FOR  ISENTROPIC  FXTERNAL  -  COMPRESSION  INLET  (SEE  REF  40) 
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FIG  75 


BOUNDARY  LAYER  BLEED  REQUIRED  TO  PREVENT  SEPARATION 
ON  SINGLE-RADIUS  TWO-DIMENSIONAL  CURVED-SURFACE 
MODEL  AT  MACH  NUMBER  2.5 

aT=  30  DEG 

FLOW  BLED  BETWEEN  as=l2  AND  25  DEG 


RADIUS  OF  CURVATURE  OF  MODEL  R 

BOUNDARY  LAYER  MOMENTUM  THICKNESS  AT  STATION  o’  80 
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FIG  76 


BOUNDARY  LAYER  BLEED  REQUIRED  TO  PREVENT  SEPARATION 
ON  SINGLE-RADIUS  TWO-DIMENSIONAL  CURVED -SURFACE 
MODEL  AT  MACH  NUMBER  3.0 

qt=  30  DEG 

FLOW  BLED  BETWEEN  a  s  =  12  AND  25  DEG 
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EFFECT  OF  WALL  TEMPERATURE  ON  REDUCED 
BOUNDARY  LAYER  SHAPE  PARAMETER 
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